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PREFACE 


The  work  reported  herein  was  performed  by  the  Sikorsky  Aircraft  Division 
of  United  Technologies  Corporation  under  Contract  DAAJ02-76-C-0066  for  the 
Applied  Technology  Laboratory,  US  Army  Research  and  Technology  Laboratories 
(AVRADCOM),  Fort  Eustis,  Virginia.  The  work  was  carried  out  under  the  tech¬ 
nical  cognizance  of  Mr.  Tom  White  of  the  Applied  Technology  Laboratory. 

The  program  was  conducted  under  the  management  of  Sikorsky  Aeromechanics 
Branch  Manager,  Mr.  Peter  Arcidiacono  and  Aeromechanics  Chief  of  Aerodyna¬ 
mics,  Mr.  J.  Rorke.  Aeromechanics  personnel  involved  directly  with  the 
analysis  and  program  development  was  Mr.  Robert  E.  Studwell. 


SUMMARY 


The  computer  program  developed  for  this  effort  is  referred  to  as  the  Heli¬ 
copter  Dynamic  Performance  (HDP)  Program.  It  was  designed  to  permit  analy¬ 
sis  of  the  capabilities  of  any  definable  single  rotor  helicopter  from  any 
heliport  location  with  specified  departure  or  approach  procedures.  Provi¬ 
sions  for  partial  or  total  engine  malfunction  at  any  point  along  the  flight 
path  are  incorporated.  The  HDP  program  generates  the  height-distance-time 
relationships  for  any  specified  vehicle  operating  weight  as  functions  of 
horsepower  input,  tip  path  plane  attitude,  and  rotor  rpm.  Provision  for 
collective  stick  input,  in  terms  of  blade  pitch  at  75%  radius  (Theta  75), 
in  lieu  of  rotor  rpm,  is  provided. 

To  accommodate  emergency  conditions,  provisions  were  made  in  the  program 
for  evaluation  of  auxiliary  power  units,  rotor  tip  rockets,  JATO,  RATO, 
fly  wheels,  drag  brakes  and  dropping  of  external  loads. 

Provisions  have  also  been  incorporated  for  calculating  the  height-veloci ty 
envelope  for  user  specified  contact  velocities  and  blade  structural  limit 
cri teria . 

The  HDP  program  is  coupled  to  the  Sikorsky  general  plot  package  to  permit 
interactive  capability  for  rapid  flight  path  evaluation.  The  program  pro¬ 
vides  an  interface  between  the  Engineering  routines  and  the  plot  package 
which  permits  the  plotting  of  40  pertinent  output  parameters.  The  plot 
interface  was  designed  so  that  any  parameter  can  be  plotted  against  any 
other  parameter  in  the  list. 

i 

In  addition  to  these  plots,  a  perspective  view  of  the  flight  path  can  also 
be  obtained.  The  user  has  the  capability  of  changing  the  viewing  position 
to  any  position  in  virtual  space.  A  plot  of  the  height-velocity  curve  for 
the  rotorcraft  being  evaluated  can  also  be  obtained  through  the  graphics 
interface  package. 

This  volume  deals  with  the  engineering  aspects  and  basic  rotor  performance 
equations  used  for  obtaining  the  force  and  energy  balance  required  during 
each  finite  time  step  of  the  overall  flight  path. 

A  parametric  trade-off  investigation  was  conducted  to  establish  the  sensi¬ 
tivity  of  various  rotorcraft  design  parameters.  This  investigation  indi¬ 
cated  that  the  parametric  sensitivity  is  highly  dependent  on  the  baseline 
configuration.  It  also  pointed  out  that  a  high  interdependence  of  basic 
design  criteria  exists.  A  means  of  evaluating  the  sensitivity  and  inter¬ 
dependence  criteria  was  established  which  can  be  used  as  a  guide  for  design 
considerations. 

A  comparative  analysis  between  program  calculated  performance  and  flight 
test  data  was  conducted  to  establish  program  accuracy. 
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HOP  PROGRAM  OVERVIEW 


The  HDP  Program  provides  a  means  of  calculating  the  space-time  relation¬ 
ships  of  a  rotorcraft  in  dynamic  (accelerating)  flight.  The  program 
basically  performs  a  step  function  integration  of  rotorcraft  unbalanced 
forces  and  energy  levels  to  establish  incremental  space  displacement. 

These  displacements  are  summed  for  a  series  of  finite  time  steps  to 
establish  the  flight  path.  A  basic  flow  chart  depicting  program  flow 
is  shown  in  Figure  1.  Typical  program  operation  proceeds  in  the  follow¬ 
ing  manner.  Given  the  design  parameters  of  the  rotorcraft  to  be 
analysed,  defined  with  23  basic  inputs  (see  Volume  2,  Users  Manual),  and 
the  initial  flight  conditions,  also  defined  by  program  input,  the  pro¬ 
gram  will  compute  the  energy  state  of  the  rotorcraft.  This  energy 
state  assumes  the  rotorcraft  to  be  in  a  steady  state  flight  condition 
at  the  initial  condition  input  values.  Additional  program  inputs 
provide  a  means  for  the  user  to  specify  variations  in  rotorcraft  velocity, 
tip  path  plane  attitude,  rotor  RPM  and  engine  power  levels.  A  series  of 
10  sequences  to  a  maximum  of  300  time  steps  is  permitted.  A  varia¬ 
tion  of  any  one  or  any  combination  of  these  inputs  will  create  an  un¬ 
balanced  force  on  the  system.  Tnis  unbalanced  force  is  computed 
assuming  the  rotorcraft  to  be  in  steady  state  flight  for  the  conditions 
which  exist  at  the  mid-point  of  each  time  step.  Since  the  rotorcraft 
acceleration  is  unknown  the  following  process  is  used  to  attain  a  solu¬ 
tion. 

Initially  the  program  assumes  the  acceleration  energy  to  be  0  and  com¬ 
putes  the  unbalanced  forces  acting  on  the  system  due  to  the  pertur¬ 
bation  on  the  input  sequence  list.  Based  on  these  unbalanced  forces 
an  acceleration  is  computed.  This  acceleration  revises  the  energy 
state  of  the  system  along  with  the  velocity  at  which  the  unbalanced 
forces  are  generated.  The  new  velocity  used  is  based  on  that  which 
exists  at  the  mid-point  of  the  time  step.  This  new  energy  state  is 
then  fed  back  and  used  to  recalculate  the  rotor  forces,  unbalanced 
forces  and  resulting  accelerations.  This  process  is  repeated  for  each 
time  step  until  the  resultant  acceleration  between  two  successive 
iterations  is  within  0.3  ft/sec  .  To  minimize  the  convergence  itera¬ 
tions  required,  each  time  step  uses  the  resultant  acceleration  frcm  the 
previous  step  as  a  starting  point. 

The  flight  path  is  obtained  by  summing  the  displacements  of  each  time 
step . 

A  discussion  of  the  equations  used  for  computing  the  rotorcraft  energy 
state,  rotor  forces,  acceleration  and  displacements  are  covered  in 
this  report. 
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ROTOR  PERFORMANCE  MODEL 


The  Helicopter  Dynamic.  Performance  Method  is  a  semi -empirical  procedure 
for  evaluating  rotor  and/or  rotorcraft  performance.  The  equa.ions  derived 
for  this  analysis  were  normalized  to  permit  performance  calculations  in 
any  flight  regime  for  any  predefined  rotor  system.  All  equations  are 
continuous  between  flight  regimes. 

The  analysis  was  based  on  an  energy  method  approach  modified  to  account 
for  blade  interference  effects.  The  effects  of  Reynolds  Number,  Mach 
Number, and  Skewed  flow  are  also  accounted  for  in  the  analysis. 


OUTLINE  OF  THEORY 


Investigation  of  several  means  of  approach  to  performance  analysis  indi¬ 
cated  that  the  energy  method  was  better  suited,  being  more  amenable  to 
analysis  as  flight  conditions  approach  zero  speed. 

The  energy  method  basically  states  that  the  power  required  for  a  rotor 
system  to  develop  a  given  thrust,  in  steady-state  flight,  is  equal  to  the 
Z  rotor  induced  power,  blade  profile  power,  airframe  drag  power,  climb 
power, and  rotor  stall  power.  The  basic  rotor  system  as  used  for  this 
analysis  is  snown  in  Figure  2. 


Figure  2.  Basic  Rotor  Force  System. 

hp  -  hp.  +  hp  +  hp  +  hp  +  hp  (1) 

a  i  o  P  c  s 
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It  can  be  seen  from  the  power  equation  that  the  use  of  the  energy  method 
reduces  itself  to  the  determination  of  each  power  component.  The  analysis 
used  for  the  determination  of  each  of  these  components  is  summarized  here¬ 
in. 


INDUCED  HORSEPOWER 


Impul 5e-Momentuin  theory  states  that  the  thrust  produced  by  an  actuator 
disc  can  be  mathematically  represented  by  the  expression: 

T  =  pnR2  ( VR  +  u)  2u  (2) 

When  VD  -  0.  u  =  [T/'2*pR 2]1/2  =  u  (3) 

K  0 

Letting  u  =  u/u„  and  V  =  Vc/u  . 

0  h  0 

Wald  esbablished  that  for  any  flight  condition 

V  =  -u  sin  aw+  (u2  sin2aw-  u2  <■  1 /u2 )  ^/2  (Reference  1)  (4) 

These  equations  are  based  on  the  following  assumptions: 

1.  The  actuator  disc  consists  of  an  infinite  number  of  blades 
uniformly  accelerating  the  air  through  the  disc  with  no  loss 
of  thrust  at  the  blade  tips. 

2.  The  power  required  to  produce  the  thrust  is  represented  only  by 
the  axial  kinetic  energy  imparted  to  the  air  composing  the  slip¬ 
stream.  A  frictionless  fluid  is  assumed  so  that  no  blade  fric¬ 
tion  or  profile  losses  exist.  Rotational  energy  imparted  to 

the  slipstream  is  ignored. 

3.  The  disc  is  infinitely  thin  so  that  no  discontinuities  in 
velocity  occur  on  the  two  sides  of  the  disc. 

In  the  actual  flow  past  a  rotor  with  a  finite  number  of  blades  there  is  a 
radial  contraction  and  a  rotational  motion  imparted  to  the  slipstream. 

This  gives  a  corresponding  increase  in  the  axial  velocity  and  a  decrease 
in  the  tangential  velocity  components  acting  on  the  blade.  If  it  is  as¬ 
sumed  that  the  "ideal  rotor"  (impulse  momentum  theory)  is  100%  efficient 
in  generating  the  induced  velocity  (u),  the  increase  in  the  induced  veloc¬ 
ity  due  to  the  contraction  of  the  slipstream  in  the  rotor  disc  plane  can 
be  related  with  the  expression: 

uactual  rotor  u i deal  rotor/B 


'Coleman,  R.  P. ,  Fpingold,  A.  M. ,  and  Stempin,  C.  W. ,  NACA  ARR  No.  L5E10. 


Where  B  represents  the  efficiency  of  the  actual  rotor  with  respect  to 
the  "ideal  rotor"  (see  section  on  development  of  the  rotor  effi 
ciency  equation) . 

Thus  it  can  be  stated  that 

u  =  [T/2ttpR2]1/2/B~ 

°actual  rotor 


or 

-  r  1  i1/2 

U.  ”  L - = - pJ 

°actual  rotor  2ttp(BR) 


(b) 


This  states  that  the  induced  velocity  generated  by  an  actual  rotor  of 
radius  R  will  be  the  same  as  that  generated  by  an  "ideal  rotor"  of  radius 
BR. 

The  horsepower  required  in  generating  this  induced  velocity  can  be  ex¬ 
pressed  as 

hip -  LL. 

1  550 


2  2 

Multiplying  by  uQ/uo  and  nondimensionalizing  thrust  with  puR  (cR)  yields 
the  coefficient  of  induced  power  as 

cp  =  550  hPi  =  Ct3/2  _  (6) 

•\ - 5 - -  —  u  A 

pirR^(oR)-5  VZ  B 

-This  equation  establishes  a  means  of  computing  the  induced  power  require¬ 
ments  for  a  single  rotor.  If  tandem  rotors  are  assumed,  further  considera 
tions  for  rotor  interference  effects  on  the  induced  power  must  be  made. 

TANDEM  AND  COAXIAL  ROTO  1  HOVER  PERFORMANCE 

Consider  two  copUnar  overlapped  rotors  as  shown  in  Figure  3. 


figure  3.  Tandem  Rotor  Overlap  Geometry. 


If  the  profile  power  is  assumed  to  be  independent  of  the  overlap  for  a 
giver,  thrust,  then  for  a  constant  induced  power  the  thrust  must  vary  in¬ 
versely  as  the  induced  velocity  for  P.  =  Tu..  Letting  T  equal  the  total 

thrust  of  both  rotors  and  u\  equal  the  'average'  induced  velocity,  the  re¬ 
lationship  of  thrust  and  induced  velocity  can  be  expressed  with  respect  to 
two  similar  but  non-overl apped  rotors  as 


where  the  subscript  o  refers  to  zero  overlap 
now  u.  =  [T/2pA]1/2 

where  A  -  2A1  +  A^  (from  Figure  3) 

and  u.  =  [To/^pA^^2 
o 

where  A  =  2nR2 
o 

1 13 

Therefore  T/Tq  =  (A/Aq)  '  for  constant  induced  power.  Similarly,  for  a 

constant  thrust  the  ratio  of  the  induced  powers  with  and  without  overlap 
is: 

Pi/Pi  =  [Aq/A]1/2  (7) 

o 

If  the  overlap  is  defined  as  d/2R  ,the  angle  y  (shown  in  Figure  2)  is  given 

by 

y  =  cos-1  (overlap) 
and  from  the  geometry  of  the  sketch 

A/A  =  1  -  1  ^1nT-c-°.s-l  (8) 

0  TT 


Considering  a  zero  overlap  condition  it  can  be  stated 
thrust  of  one  rotor  equal  to  half  the  total  then2p. 

’l 


that 

=  P. 

i  , 


by  assuming  the 


Where  subscript  1  refers  to  one  rotor  and  2  to  two  rotors, 
overlapped  condition 


A 

■  2Pi  £/r  > 

1 1  atot 


1/2 


Thus  for  an 

(9) 
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where  P.  is  based  on  the  thrust  and  area  of  one  rotor. 

n 

Defining  df^  as  the  total  rotor  interference  factor  in  hover  due  to  over¬ 
lap  then  the  total  induced  power  for  hover  can  be  expressed  as 


P.  =  2P.  (1  +  df jj/ 2 )  and  consequently 

7  'l  M 


Ti 


solving  for  dfH  yields 


dfu  -  2  [- 


1 


1  - 


7-] 


1/2 


-  sin-, -cosy 

7T 


(10) 


This  analysis  was  based  on  the  rotors  being  coplanar  and  thus  assumes 
that  the  stream  tube  contraction  for  the  two  rotors  is  equal.  When  the 
rotors  do  not  lie  in  the  same  horizontal  plane  the  effects  of  differential 
slipstream  contraction  should  be  considered.  Full  scale  test  data  indi¬ 
cated  that  this  effect  can  be  approximated,  for  hover  at  least,  by  multi¬ 
plying  the  coplanar  interference  factor  by  the  expression  (1+19  vs/R). 
This  term  is  strictly  empirical  and  justification  for  its  use  is  only  with 
its  apparent  agreement  with  available  test  data.  To  maintain  conformity 
with  the  actuator  disc  theory,  the  maximum  value  of  the  interference  fac¬ 
tor  (the  product  of  the  vertical  separation  and  overlap  expressions)  is 
limited  to  2.0  where  VS  refers  to  the  vertical  displacement  between  rotors. 


TANDEM  AND  COAXIAL  ROTOR  FORWARD  FLIGHT  PERFORMANCE 

A  comparison  can  be  made  between  single  rotor  and  tandem  rotor  helicopters 
at  the  same  gross  weight  by  assuming  that  the  thrust  is  equally  divided 
between  the  two  rotors  of  the  tandem.  According  to  simple  momentum  rela¬ 
tionships,  the  induced  power  in  forward  flight  for  an  isolated  rotor  is 

Pi  =  T2/2PAu 


for  the  rear  rotor  of  a  longitudinal  tandem  this  becomes 
P.  =  (T2/2PAu)  (1  +  df  ) 

'r  hr 
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and  for  the  front  rotor 

Pi(r  =  (T2/2pAu)  (1  +  dfRp) 

where  subscripts:  FR  refers  to  the  interference  of  the  front  rotor  on  the 

rear  rotor 

RF  refers  to  the  interference  of  the  rear  rotor  on  the 
front  rotor. • 


Thus  the  total  induced  power  would  be 
T2 


P.  +  P.  = 

ic  2pAu 


(2  +  df  +  df  )  or  P. 

RF  ’TOTAL 


=  2P.  1  + 

pp  nr  1  tat  «i  1  * 


dfFR  +  dfRF 


'F  FR  RF  'TOTAL  '1  2 

Reference  2  derives  an  expression  for  coplanar  longitudinal  tandems  where 


>2COS2f 

,2_2. 


df  is 

defined  as 

dfFR  =  ( 

1  + 

P2  + 

Psinq)/  1  +  P 

and 

0 

dtRF  =  ( 

1  + 

P2  - 

psi n;)/  1  +  p 

where 

C  =  TAN'1 

2pAf 

:u2/l, 

,5Tf 

The  factor  1.5  accounts  for  partial  wake  contraction  at  the  rear  rotor. 
Summing  the  terms  yields 

dfLF  =  dfFR  +  dfRF  =  2/(1  +  pZcoi2^  0!) 


To  account  for  variations  in  vertical  displacement  between  rotors  test 
data  indicates  that  this  term  should  be  adjusted  by  the  term  (1  -  1.35 
dv/R).  This  assumes  the  rear  rotor  to  be  higher  than  the  forward  rotor. 
Where  dv  is  the  vertical  displacement  between  rotors. 

While  this  expression  appears  reasonably  valid  in  high-speed  flight, its 
usage  becomes  questionable  iri  the  low  speed  regime.  To  relate  the  df  in 
forward  flight  to  that  derived  for  hover  it  was  assumed  that  a  linear 
transition  between  the  hover  point  arid  the  level  flight  curve  exists.  On 
this  basis  the  df  in  the  low  speed  regime  could  be  expressed  as: 


df 


dfH  + 


dfLF  ~  dfH 
CTAN 


(12) 


JMcCormick,  B.  W. ,  Jr.,  AERODYNAMICS  OF  V/STOL  FLIGHT,  Academic  Cress, 
1967. 
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where  is  the  wake  skew  angle  at  which  a  line  projected  from  df^  is 

tangent  to  the  df  curve  for  the  rotor  system  in  question  (see  figure  4). 


ASSUMED 

transition 


; 

Figure  4.  Tandem  Rotor  Interference  (Typical). 


The  use  of  these  equations  yields  excellent  agreement  with  test  data  for 
tandem  rotor  helicopters  with  shaft  offset  but  tend  to  breakdown  as  the 
shaft  offset  distance  (d)  approaches  zero,  (see  Figure  51 

Considering  a  two-rotor  system  with  substantial  shaft  it  can  be  seen 

that  as  forward  speed  is  increased  the  stream  tubes  ac  ••■.ed  by  each 
rotor  interfere  by  an  increasing  amount. 


Figure  5.  Tandem  Rotor  Inflow  (Typical). 
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Thus,  as  stated  by  the  equations  above, the  interference  sn  hover  is  low 
and  as  V  tends  to  infinity  df^p  approaches  2,  lessening  the  effective 
vertical  displacement. 

Now,  as  the  shaft  offset  distance  is  decreased  this  effect  tends  to  re¬ 
verse,  i .e. ,  the  interference  is  highest  in  low-speed  flight  and  decreases 
as  speed  increases  (figure  6). 


Figure  6.  Co-axial  Rotur  Inflow  (Typical). 


To  properly  account  for  the  condition  when  the  offset  is  zero  (coaxial 
rotors)  it  is  suggested  that  the  interference  effects  be  expressed  as 

df  =  dfH  cos2;  (13) 


and  defining  the  induced  power  of  each  *'otor  as 


Pi  =  [l.lT2/2pf  )  ('  +  df  ) 
upper  u 

and  P.  =  [1.4T2/2dAu]  (1  +  df  ), 

'lower  T1 

2  d1u  +  dfl 

thus  P.  =  2.5TV2pAu  [1  +  - LJ  (14) 

1  total 


These  expressions  assume  that  each  rotor  is  carrying  50%  of  the  lift.  The 
factors  1.1  and  1.4  are  purely  empirical  and  are  the  estimated  entrainment 
effects  of  one  rotor  on  the  other. 

It  must  be  stressed  again  that  interference  expressions  defined  above  are 
semi -empi rical  in  nature  and  exact  distributions  would  require  a  more 
rigorous  form  of  computation. 


PROFILE  HORSEPOWER 


The  impul se -momentum  concept  assumes  the  rotor  to  be  operating  in  a  fric¬ 
tionless  fluid  and  profile  drag  losses  are  not  accounted  for.  To  account 
for  the  profile  drag  losses,  the  following  relationships  were  established. 
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-  -  _  -  a 


The  thrust  produced  by  a  rotor  system  can  be  expressed  on  a  strip  analysis 
basis  as 


T  =  --  f2l]  dv  /  R  p/2  CU2  C,  dr 
2tt  o  o  L(M) 


Accounting  for  the  influence  of  the  reversed  velocity  region  on  the  re¬ 
treating  blade  the  integral  becomes 


T  = 


2n 


dip 


rR 


p/2  CU  C. 


•(M) 


dr  -b/ 2 7i 


.2 tt  ,  .-uRsin^  9 
,  d"  '  P/2CU2C,  dr 
0  L(M) 


This  expression  for  thrust  in  terms  of  does  not  account  for  slip¬ 

stream  rotation  and  is  therefore  representati ve  of  an  "ideal  rotor".  Uti¬ 
lizing  the  analogy  that  an  actual  rotor  of  radius  R  is  equivalent  to  an 
"ideal  rotor"  of  radius  BR,  the  thrust  developed  by  an  actual  rotor  can  be 
computed  by  integrating  along  the  blade  to  radius  BR. 


Substituting  C-j-  piiR2(nR)2  for  T 


and 


fir  +  usiR  sin>j'  for  U 


Integrating  and  solving  for  C 


KM) 


yields 


=  2  [; 


1 


'\  *"  '  /  up  L  o  ___  O 

L(M)  e  Bu -  4u 

3  2  9h 


3- 


(15) 


(Definition  of  a  is  discussed  later), 
e 

The  required  mean  lift  coefficient  C|_^  can  be  established  in  terms  of 
Cj  and  the  resulting  mean  profile  drag  coefficient,  6,  can  be  established 
from  the  6-C.  relationship  of  the  airfoil  in  question. 

l(m) 

In  Glauert's  theory  of  energy  loss  the  expression  for  profile  drag  with 
skewed  flow  taken  into  account  is 


e  6 


2, 


Do/L  “  aTET  <’  > 


Assuming  L  =  T  cos  a  and  cos  a  =  1 


19 


I 


and  substituting  T/pttR2(oR)2  for  CT 
and  V/f;R  for  u 

yields  D QV  =  (0  6pnR^  (fiR)3/8l  (1  +  nu2) 

expressing  in  nondimensional  terms  yields 

°e5  ? 

Cp  =  -p-  (1  +  nu  ) 


or  in  terms  of  horsepower 


hP0  = 


°e5 


(1  +  n|j  ) 


2  \  p  7t  R  2  ( >  R ) 2 


550 


(16) 


07) 


Where  n  is  a  function  of  p  determined  by  the  equation: 

1  +  np2  /  '  (x2  +  2xpsin^+p2)3/2dxd. 

Where  =  V/ uR  <  1  ,  x  =  r/R 

Lacking  a  simple  integral  of  this  equation,  Glauert's  method  (Reference  3) 
of  obtaining  approximations  by  evaluating  the  inner  integrals  at  the  points 
of  ^  =  0,  90°,  180°  and  270°  and  averaging  has  been  adopted.  Using  this 
procedure,  numerical  values  of  u  and  n  were  evaluated  and  are  given  in  the 
following  table. 


1* 

0 

— 

.2 

n 

.5 

.6 

1  .0 

n 

4.50 

n 

4.63 

4.73 

BE9 

5.03 

5.22 

6.13 

The  mean  profile  drag  coefficient  (-S)  is  related  to  the  mean  lift  coeffi¬ 
cient  (C.  )■  This  relationship  is  generally  related  as  a  cubic  function  of 

(M) 

C,  .  A  typical  relationship  for  an  NACA  0012  airfoil  is  shown  in  Figure  7. 

(M) 

DEVELOPMENT  OF  THE  ROTOR  EFFICIENCY  EQUATION 

The  previous  sections  show  that  the  calculations  of  the  induced  and  profile 
powers  are  dependent  on  establishing  th  otor  efficiency.  The  following 
analysis  is  used  to  establish  the  rotor  efficiency. 

’Glauert,  H. ,  A  GENERAL  THEORY  OF  THE  AUTOGYRO,  ARC  R&M  No.  Ill,  1926. 
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Simple  momentum  theory  states  that  the  thrust  developed  from  an  actuator 
disc,  operating  in  the  normal  working  state,  can  be  determined  from 
equation  2. 

In  the  actual  flow  past  a  rotor  there  is  a  radial  contraction  and  a  rota¬ 
tional  motion  imparted  to  the  slipstream.  This  gives  a  correspondi ng  in¬ 
crease  in  the  axial  velocity  and  a  decrease  in  the  tangential  velocity 
component  acting  on  the  blade.  Representing  a  hovering  rotor  with  an  in¬ 
finite  series  of  airfoil  sections  operating  in  one  plane,  the  relative 
velocities  acting  on  each  airfoil  section  can  be  represented  as  shown  in 
Figure  8.  The  distance  between  consecutive  airfoil  sections  would  be 
2rr/b  the  induced  velocity  would  be  u  and  the  whole  system  would  be  moving 
at  a  velocity  ,Qr. 


Figure  d.  Blade  Interference  Effects. 


Let  us  consider  a'  as  the  axial  interference  factor  and  a"  as  the  rotation 
al  interference  rotor.  Then  u  would  become  u(l  +  a')  and  ;>r  would  be 
3r(l-a").  The  forces  developed  by  the  airfoil  section  due  to  these  veloci 
ties  can  be  represented  as  shown  in  Figure  9.  The  solid  lines  represent 
the  "Ideal"  rotor  while  the  dash  lines  represent  the  'Actual ’  rotor. 


‘  i  : 

| 

i — -  - - 1 


Figure  9.  Interference  Effects  on  Blade  Forces. 
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It  can  be  seen  from  Figure  9  that  the  vertical  and  horizontal  forces,  Cy 
and  Cx,for  the  Ideal  rotor  element  can  be  determined  from  the  equations 

0^  =  coso  -  CQ  sin-; 


Cx  =  CL  sini+Cp  cosy 


4>  -  tan-1  u/ar 

And  for  the  Actual  rotor  element 


C  '  =  C.  '  cos^; '  -  Cn'  sin^.' 
y  l  0 

Cx'  =  CL‘  sinf'  +  Cp  cos y 

<?■’  =  tan'1  u/  f:r  ( ( 1  +  a‘ )/ (1  -a"  ) ) 

If  the  actual  blade  element  is  to  maintain  the  same  lift  coefficient  as  the 
ideal  blade  element  the  a'  must  equal  a.  If  a'  is  to  equal  a,  6  must  be 
increased  by  the  amount  $'  -  When  the  ideal  and  actual  blade  elements 

are  operating  at  the  same  angle  of  attack,  they  will  both  develop  the  same 
lift  and  drag  coefficients  but  the  vertical  and  horizontal  components  will 
not  be  equal  as  ■*>’  >  <j>.  In  order  for  the  actual  rotor  to  develop  the  same 
Cy  and  Cx  values  as  the  ideal  rotor  the  angle  of  attack  must  be  increased 
so  that  a 1 >  a.  Since  the  actual  rotor  must  maintain  a  higher  angle  of 
attack  to  attain  the  same  vertical  force,  a  higher  CL  value  must  be  de¬ 
veloped  which  will  give  a  corresponding  increase  in  C^. 

From  simple  momentum  theory  the  induced  velocity  (u)  can  be  expressed 
from  equation  2  as 


u 


~V/Z 

2pnR 


(1«) 


Thus  when  V  =  0  the  induced  velocity  thru  an  ideal  rotor  system  can  be 
expressed  as 


u 


c 


] 


1/2 


(3) 


23 


Dividing  equation  18  by  equation  3  yields 


u/u0  =  -  art  n/4  (^)2  + 1]1/2  =  u  (i9) 

0  0 

Thus  when  V  =  0,  u  =  1  for  an  Ideal  rotor. 

In  Figure  8  let's  assume  that  the  first  blade  element  represents  the  ideal 
rotor  and  the  second  element  an  actual  rotor.  On  this  basis  it  can  be 
seen  that  the  ratio  of  induced  velocity  of  an  actual  rotor  to  an  ideal 
rotor  would  be  greater  than  1,  since  u(l  +  a')>uo.  This  indicates  that  an 

actual  rotor  is  less  efficient  than  an  ideal  rotor. 

The  efficiency  of  any  real  rotor  system  can  be  evaluated  from  the  relation¬ 
ship 

u(1  +  u')B  =  u  =  (20) 

0  2pnR* 

Where  B  represents  the  rotor  efficiency. 

This  equation  can  be  rewritten  as 

u(l  +  a')  -  [ - (21) 

2p:;(BRr 

Which  states  that  the  induced  velocity  through  an  actual  rotor  of  Radius  R 
can  be  represented  by  computing  the  induced  velocity  of  an  ideal  rotor  of 
radius  BR.  Based  on  this  relationship  the  rotor  efficiency  was  evaluated 
from  test  data  in  the  following  manner. 

The  induced  velocity  through  an  actual  rotor  disc  was  determined  from  the 
equation 


u( 1  +  a 1 )  =  550  hp.  /T 


(22) 


Where  hpi  -•  hpa  -  hpQ  -  hpp 


hp. 


is  the  induced  power 


hn  is  the  total  power  measured  at  the  rotor 
Ka 
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hp 


o 


is  the  profile  power 


hpp  is  the  parasite  power  required  by  the  airframe  for  wind 
(equivalent  forward  speed). 

The  profile  power  was  determined  from  the  equation 


hP, 


°e6 


(1  + 


n,2) 


pttR^o 


550 


(17) 


where  6  is  a  function  of  C. 

L(n) 

and  C.  =  2k  0,/.- 

L(M)  t  'e 

1 

where  k  =  B2  Bu2  4u2 
3  2  “  9  w 

See  previous  section  for  the  derivation  of  these  equations. 


The  use  of  the  rotor  efficiency  B  in  equation  15  accounts  for  the  in¬ 
creased  lift  (C^  )  requirements  and  the  corresponding  increase  in  profile 

drag  (5)  for  the  actual  rotor. 


Total  rotor  losses  B  were  then  determined  through  an_iterative  solution  of 
these  equations.  The  initial  assumption  being  that  B  =1  in  equation  (15) 
and  a  value  of  B  determined  from  equation  21.  Jlris  value  of  B  was  then 
substituted  in  equation  15  and  the  process  repeated  until  convergence. 

Utilizing  flight  test  data  obtained  from  hovering  flight,  values  of  B  were 
determined  from  equation  21  for  various  aircraft  over  a  range  of  thrusts. 
The  mean  profile  drag  coefficient  vs.  mean  lift  coefficient  polar  for  the 
NACA  0012  airfoil  was  assumed  to  be  substantially  correct  and  used  in  this 
analysis.  The  calculated  values  of  B  showed  variations  depending  on  the 
number  of  blades,  blade  twist  and  thrust  coefficient. 


The  values  of  B  calculated  for  a  five  bladed  rotor  with  a  blade  twist  of 
-14°  were  approximated  by  the  equation 

[1  ,34Ct]1/5 

D  -  1  * 


This  was  then  written  in  a  general  form  as 


B  = 


[1 ,34CT]1/b 


b 
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(23) 


It  should  be  noted  that  this  equation  form  yields  B  -  1  as  b  ->■  **■  which 
must  be  the  case  as  the  rotor  system  approaches  the  actuator  disc  form  of 
the  ideal  rotor.  Applying  equation  23  to  rotors  with  three  and  four 
blades  showed  that  the  test  data  were  displaced  from  the  curve  predicted 
by  equation  23.  It  was  also  noted  that  this  displacement  was  constant  for 
a  given  blade  twist.  This  displacement  used  in  conjunction  with  the  cor¬ 
responding  blade  twist  was  approximated  with  the  equation 

aB,  =  -  [0. 143250  +  0.035] 

D  ^  C 


where  e  is  the  blade  twist  expressed  algebraically  ir,  radians. 


A  general  equation  for  the  total  rotor  losses  in  a  hovering  condition  was 
then  written  as 


[1.34Ct] 


1/b 

-  [0. 14325a  + 

e 


0.035] 


(24) 


In  translational  flight  the  value  of  V  is  not  zero  and  the  va|ue  of  u  is 
not  one.  Walds  equation.  Reference  1,  derives  a  function  of  V  vs.  u  for 
various  rotor  disc  angles.  In  translational  flight  an  ideal  rotor  should 
yield  values  of  u  obtained  from  this  equation  atthe  proper  rotor  disc 
angle  and  V  value.  Equation  4  yields  vaUies  of  V  that  would  be  obtained 
on  an  ideal  rotor.  To  obtain  values  of  V  that  are  obtained  on  an  actual 
rotor, equation  19  must  be  rewritten  as 


V/u 


o 


v/C — 

2Ttp(BRr 


1/2 


The  total  rotor  losses  in  translational  flight  can  then  be  determined 
through  an  iterative  solution  using  equations  21,  22,  17,  15  and  25. 


(25) 


Flight,  test  data  were  used  in  determining  values  of  B  in  translational 
flight.  The  calculated  values  of  B  showed  variations  with  the  number  of 
blades,  blade  twist,  thrust  coefficient  and  translational  speed.  When  the 
B  value  as  determined  from  equation  24,  hovering  flight,  was  subtracted 
from  the  B  value  determined  from  translational  flight  tests  it  was  noted 
that  the  difference,  'A Bp ,  was  constant  at  a  given  CT  and  »  value  for  all 
aircraft.  It  was  also  noted  that  aBu  varied  with  CT  and  u- 


The  coning  angle  of  the  rotor  blades  increases  as  the  weight  of  the  air¬ 
craft  increases.  This  coning  angle  is  constant  for  any  given  weight, 
temperature,  altitude  and  rpm.  Flight  tests  indicate  that  for  a  given  set 
of  conditions  the  coning  angle  decreased  as  the  translational  speed  in¬ 
creased,  up  to  the  minimum  power  speed,  and  then  increased  with  a  further 
increase  in  speed.  This  decrease  and  increase  of  coning  angle  in  transla¬ 
tional  flight  can  be  attributed  to  changes  in  blade  loading,  as  shown  in 
Figure  10. 
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figure  10.  Effect  of  Thrust  Distribution  on  Blade  Coning. 
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These  variations  in  blade  load  distribution  indicated  that  the  losses  in 

the  system  should  also  vary  with  speed.  Noting  that  the  term  V/u  can  be 

1/2  0 
expressed  as  u(2/CT)  and  that  the  blade  loading  distribution  will  also 

vary  with  azimuth  position,  a  generalized  expression  for  the  variation  in 
rotor  efficiency  with  translational  speed  was  written  as 

2  ,  _  2 


ABu  =  xw(2/Ct),/2  +  ip-  +  ij1/2  -  ip-  +  1] 


1/2 


(26) 


4,  0°  &  180°  90°  270° 

where  x,  v  and  z  are  unknown 

Values  of  aBu  were  determined  by  subtracting  the  B  values  calculated  from 
equation  24  from  the  B  values  calculated  from  translational  flight  at 
various  test  values  of  u  and  C^.  Calculating  for  several  conditions  and 

substituting  the  results  in  equation  26  yielded  values  of  x,  y  and  z  as 
x  =  0.0905,  y  =  0.5,  z  =  0.6974 
Equation  14  then  becomes 


iB^  =  0.0905u  Cl— 31 /2 
^T 


1/2 


-  L 


0 . 6974  l 


+  1] 


1/2 


(26a) 


27 


The  total  expression  for  the  rotor  losses  due  to  slipstream  contraction 
slipstream  rotation  and  blade  tip  losses  in  any  flight  condition  was  then 
written  as 


B  = 


[1.34CT]1/b 

b 


+  0. 0905m  lU]/Z  +  LjJr-  + 

(27) 

2 

.  [Mi74u_  +  l]1/2  -  [0.14325  +  0.035] 

V  J  V 


The  equation  for  B  having  been  developed  by  comparing  the  induced  velocity 
through  an  actual  rotor  to  the  induced  velocity  through  an  ideal  rotor  is 
then  in  effect  the  efficiency  of  the  actual  rotor.  The  increased  pro¬ 
file  drag  requirements  due  to  the  higher  lif„  needed  to  maintain  the  same 
Cy  force  as  the  ideal  rotor  are. a.lso  accounted  for  in  this  analysis. 

.  ««  fr¬ 

it  should  be  noted  that  in  the  development  of  B  the  flight  path  velocity 
was  used  in  determining  values  of  p.  On  this  basis  the  resultant  flight 
path  velocity  must  be  used  for  the  determination  of  p  in  the  rotor  effi¬ 
ciency  equation,  i .e.  , 


PB  =  [VH2  +  vv2]1/2/siR 


(28) 


THE  INFLOW  RATIO  IN  DESCENT 


The  power  requirements  of  a  rotor  system  in  low-speed  flight  are  largely 
dependent  on  the  induced  power  or  more  specifically  the  induced  velocity 
through  the  rotor.  In  order  to  fully  utilize  the  equations  derived  a 
means  of  determining  the  induced  velocity  during  descent  must  be  estab¬ 
lished.  This  is  necessitated  as  'Walds’  equation,  derived  in  Reference  1 
is  deemed  not  applicable  in  the  negative  velocity  range  between  the  values 
of  V/uQ  =  0  and  V/uQ  =  -2. 

Investigations  into  the  vertical  descent  range  V/uQ  =  0,  -2  have  been  con¬ 
ducted  by  Hafner,  Glauert  and  Lock  who  have  proposed  various  transition 
curves  across  this  range.  Various  presentations  proposed  are  summarized  in 
Reference  4.  Reference  5  also  proposes  a  transition  curv.j  based  on  test 
data  obtained  from  Reference  6.  It  should  be  noted  that  LheU/u0,  V/uQ  re¬ 
lationships  for  the  proposed  transition  curves  are  based  on  "ideal"  rotor 

“Shapiro,  J.  ,  PRINCIPLES  OF  HELICOPTER  ENGINEERING,  McGraw  Hill  Book 
Company,  Inc.  ,  1955. 


■Gessow,  A.,  and  Myers,  G.  C.,  Jr.,  AERODYNAMICS  OF  THE  HELICOPTER, 
MacMi Ilian  Company,  1952. 

r,NACA  MR  No.  L5U09a. 
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concepts,  i.e.,  no  losses. 


On  this  basis 


u/uo  = 
where  hp 

a 

hp 

o 


550(hpa  -  hp  ) 

- and  V/u  =  V./LT/2.PR2]17 

T  [T/2npR2l1/2  0  R 

measured  from  flight  test 

computed  for  flight  test  conditions 


c 


T  =  GW 

V^  Resultant  Flight  Path  Velocity 

The  analysis  conducted  previously  modifies  these  equations  to  account  for 
rotor  losses  and  states  that  for  actual  rotors  the  equations  would  be 


u/u 


o 


550Cipa  -  hpQ  ) 

T  [T/2q*(BR)2]1/2 


and  V/u 

o 


[ - b-il 

2pir(BR; 


4 


1/2 


where  T  =  f (HP) ;  hpo  =  f(B) 


Reducing  the  test  data  for  near  vertical  descent,  as  given  in  Reference  6  on 

this  basis,  yielded  u/uQ  =  1.0  for  all  the  test  points*  This  analysis 

used  the  rotor  efficiency  equation  for  the  B  term  in  these  equations.  These 
results  indicate  that  the  transition  curve  between  V/uQ  =  0  and  V/uQ  -  -2 

can  be  expressed  on  an  actual  rotor  basis,  by  u/uQ  =  1.0.  The  results  of 
this  analysis  are  shown  in  Figure  11. 

Although  some  scatter  is  apparent  in  the  results,  it  is  interesting  to  note 
that  the  maximum  variation,  when  converted  to  power,  is  only  10  horsepower 
for  the  data  plotted  as  squares  in  Figure  11.  The  maximum  variation  for  the 
hover  data,  plotted  as  circles  is  3.5  horsepower. 

These  results  then  define  the  valid  boundaries  of  'Walds'  equation.  These 
boundaries  would  be  between  a  =  0°  and  =  +90°  for  all  positive  values  of 
V/uo  including  V/uQ  =  0  and  for  a  =  -90°  for  all  V/u  values  less  than  -2. 

For  the  flight  regime  between  V/uQ  =  0  to  -2,  u/uQ  =  1.  Utilizing  these 

considerations  the  relationship  of  u/uQ  to  V/u  ,  as  determined  from  'Walds' 

equation,  within  the  specified  limitations,  is  plotted  in  Figure  11.  These 
boundaries  are  used  for  this  analysis. 


*  The  four  test  points  shown  in  Reference  5  for  V/u  values  between  -0.1 
and  -0.4  could  not  be  found  in  Reference  6. 


RATIO 


Figure  11.  Induced  Velocity  Ratio  as  a  Function  of  Trans lational 
Velocity  Ratio. 


While  these  boundaries  establish  a  means  of  computing  the  induced  velocity 
during  steep  descents,  it  does  not  establish  the  maximum  rate  of  descent 
for  which  controllability  can  be  maintained.  Assuming  that  control  can  be 
maintained  up  to  the  onset  of  the  vortex  ring  state,  then  the  limit  condi¬ 
tions  can  be  determined  by  defining  the  vortex  ring  boundaries.  The  defi¬ 
nition  of  this  boundary  is  based  on  the  following  analogy. 

Consider  a  rotor  in  a  hovering  condition  with  the  rotor  producing  a  uniform 
induced  velocity  U  .  The  stream  tube,  produced  by  the  rotor,  is  surrounded 

by  a  tube  of  vorticity  generated  from  the  blade  tips.  For  this  condition, 
whirl  stand  measurements  have  indicated  that  the  rate  at  which  the  centers 
of  the  vortex  cores  move  downward  is  approximately  52.5%  of  the  average 
induced  velocity  U  .  Assuming  that  the  onset  of  the  vortex  ring  occurs 
when  the  rate  of  descent  of  the  rotor  increases  to  the  point  where  the 
vortex  cores  no  longer  move  away  from  the  rotor  disc,  then  the  critical 
rate  of  descent  would  be 


ROD 

CRIT. 


-.525  u 

o 


In  order  to  substantiate  the  assumptions  of  this  analysis,  flignt  tests 
were  conducted  at  Sikorsky  Aircraft,  with  an  S-61N  rotorcraft.  These 
tests  were  set  up  to  explore  this  flight  regime  and  to  determine  the  maxi¬ 
mum  velocity  ratio  for  which  positive  rotor  control  can  be  maintained. 

The  maximum  velocity  ratio  tested,  which  resulted  in  an  extremely  hard 
landing  and  landing  gear  damage  (V/u  =  -.68),  is  shown  in  Figure  11.  One 
other  test  point  which  resulted  in  amoderately  hard  landing  with  no  rotor- 
craft  damaqe,  V/uQ  =  -.528,  is  also  included  in  this  figure.  In  all, 

623  near  vertical  descents  were  conducted  in  the  V/Up  regime  ranging 
from  -.05  to  -.68.  However,  only  the  two  extreme  points  were  considered 
pertinent  for  the  purpose  of  this  section. 

The  results  of  these  flight  tests  indicated  that  positive  rotor  control 
could  be  safely  maintained,  provided  the  rate  of  descent  did  not  exceed 
50%  of  the  computed  induced  velocity.  On  this  basis  it  was  stated  that 
operation  within  the  "Vortex  Ring"  region  of  flight  (V/uQ  =  0  to  -2)  can 

be  safely  conducted  in  close  proximity  to  the  ground,  provided  that  the 

maximum  rate  of  descent  at  any  point  along  the  flight  pai.  dots  not  exceed 

.5u  .  These  results  indicate  that,  while  the  rotor  airflc  is  assumed  not 
o 

to  breakdown  until  V/uq  =  -.525,  positive  rotor  control  can  only  be  main¬ 
tained  up  to  a  V/uQ  -  -.5.  The  positive  rotor  control  limit  (V/u  Q  _  ^ 

is  used  for  determining  operating  capabilities  in  the  Helicopter  Dynamic 
Performance  Program. 


31 


BLADE  TWIST 


The  blade  twist  basically  establishes  the  spanwise  load  distribution  which 
in  turn  affects  the  efficiency  of  the  rotor  system.  The  analysis  conducted 
herein  has  been  based  on  linear  twisted  blades.  To  approximate  the  effects 
of  nonlinear  twist  on  rotor  efficiency  the  twist  distribution  is  converted 
to  a  linear  twist  equivalent  in  a  manner  similar  to  establishing  the  effec¬ 
tive  solidity.  The  expression  used  for  estimating  the  equivalent  linear 

twist  0  is 
e 


ee(x-0.8) 


x3d  x  ; 


6  X 


x3d  x 


(29) 


where  x  -  r/R 

“  denotes  reference 


to  0  twist  at  x  =  0.8 


PARASITE  HORSEPOWER 


To  determine  the  power  required  to  overcome  the  airframe  drag,  the  total 
drag  force  was  divided  into  horizontal  and  vertical  components. 

Initially  dealing  with  the  horizontal  component  and  expressing  drag  in 
general  terms  as 


n  _  f p» 

0  -  ~y~ 


(30) 


2  7 

And  nondimensionalizing  with  the  thrust  factor  pnR  (rtR)  yields 


D/PttR2(p.R)2 


FpV2  =  _F_2p2  =  c 

2pttR2(pR)2  2jiR  DP 


Expressing  the  nondimensional  parasite  drag  in  power  coefficient  terms 
yields 


or  in  terms  of  horsepower 

hp  r  p  nR2 ( uR ) 3 

p  "  LPp  550 


(31) 


(32) 
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VERTICAL  DRAG 


The  vertical  component  of  drag  is  generally  added  to  the  operating  weight 
to  determine  the  rotor  thrust  required.  (Steady-state  flight.)  Utilizing 
the  general  expression  for  drag, 

°v  ‘  CD  S= 

7T 

Determining  in  terms  of  rotor  disc  area 


arid  V  as  u  +  Vy 

2 

Then,Dy  =  CQ  (U  +  Vy)2 


Letting  CQ  -  / 2 

V  ir 

2  2 

and  nondimensional i zi ng  with  pnR  (uR)  yields 
□y/PT,R2(siR)2  -  C  (u  +  Vy)2/(^R)2 

Multiplying  by  uQ/uo  and  appropriately  expressing  uQ  as  Cj  u  (.qR)/v/2  B, 
the  expression  can  be  nondimensional i zed  to 

o  ?  2Bu  9  CT 

D  /pnR^fiR)^  =  Cn  [u/u  +  ~^rppT  ~Zf 
v  Uy  0  CT  '  ^  2B*" 


The  expression  for  the  vertical  component  of  thrust  in  nondimensional 
terms  then  becomes 


Cj  COS  3jpp 


-  Cw  +  CQv[u/u0  + 


/?■  B 


-.2  W_ 

1/2J  2g2 


(33) 


When  canted  tail  rotors  are  considered  the  vertical  component  of  the  tail 
rotor  thrust  must  be  subtracted  from  this  equation. 

The  term  Cn  can  be  defined  as  the  vertical  drag  coefficient. 

UV 


9 
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Defining  the  vertical  drag  as  a  percentage  of  operating  weight  with  the 
helicopter  in  an  out-of-ground  effect  hover,  the  term  CQ  can  be  defined 

from  the  equation 

°v  '  \  ‘A 

2 

Using  tt  R  for  S 

3  TT 

and  uQ  for  V 

Substituting  equation  (5)  for  uQ  and  letting  GW  +  NGW  =  T 

where  N  is  equal  to  the  percentage  of  operating  weight  for  vertical  drag 
then 


N  GW  =  Cn  p 
UV 


GW+  NGW  d2 
3-  ™ 

2roB  R 


Solving  for  Cn  yields 
UV 

\  *  2“2  ('  -  TTN> 


(33a) 


CLIMB  POWER 

The  rate  of  climb  (ROC)  for  any  flight  condition  can  be  expressed  in  terms 
of  the  excess  power  from  the  equation: 

ROC(FPS)  -  (HPavailabu  -  HPREQUiRED)!550/T  cos»TPP)No 

Where  N  is  the  rotor  aerodynamic  efficiency  in  climb. 

0 

The  velocity  and  force  vectors  in  a  steady-state  climb  (Figure  12)  indicate 
that  the  power  required  to  sustain  a  climb  rate  (V  )  must  be  equal  to 

(T  cosr,Tpp)yv  ,  since  T  coSaTpp  =  D^siny. 

The  term  T  (rotor  thrust)  is  used  here  as  the  thrust  required  to  sustain 
the  climb  condition,  and  Vy  =  Vsiny  then  (TcosaTpp)  Vy  =  DcV  or  the  climb 

power  can  be  expressed  in  terms  of  an  equivalent  drag  force:  Climb  power 
required  =  DcV. 
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Figure  12.  Velocity  and  Force  Vectors  in  Climbing  Flight. 

Nondimensionally  the  additional  power  required  to  climb  can  be  expressed 
as 

ACq/c  =  (“CD/o^CLIMB^R//rV 

or  Nq  --  (aC^/o/ jCg/o)up 

It  should  be  noted  that  the  Nondimens ional  Velocity  (uR)  is  a  function  of 
flight  path  velocity, thus 

uR  =  [Vh2  +  Vv2]1/2/nR 


The  value  of  Nq  can  be  approximatea  by  determining  the  aCq/o  required  by 
the  aCd/c  =  f (T ,  —)  at  the  proper  u,  M,  gQ  and  «e.  Evaluation  of  the  rotor 
aerodynamic  efficiency  (No)  indicated  that  it  is  dependent  on  atmospheric 
conditions  and  its  value  can  be  reasonably  approximated  by  substituting 
the  density  ratio  (d/pq)  for  Nn.  On  this  basis  the  climb  power  require¬ 
ments  are  expressed  for  this  analysis  as 

hpc  =  (T  cos  ajpp^Vv  (34) 

2  3 

Multiplying  by  550/ p ir R  (fiR)  yields 

Up  -  (By  eos  Uyp(j)u  (p  /p)  (35) 

CLIMB  ‘  TP-  v  o 
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ROTOR  SOLIDITY 


The  solidity  of  a  rotor  having  rectangular  blades  may  generally  be  ex¬ 
pressed  as  bC/nR.  Since  this  method  of  analysis  deals  with  power  or 
torque,  the  equivalent  chord  used  for  determining  the  solidity  should  be 
based  on  torque  to  obtain  a  more  accurate  determination.  This  can  be  ex¬ 
pressed  for  any  blade  shape  and  account  for  various  root  cutout  areas  as 

C  =  /  1  Cx3dx/  /  1  x3dx  (36) 

O  0 

Substituting  into  the  general  solidity  equation  bC/nR  yields 

o  =  ~p"  T  Cx3dx/  /  1  x3dx  =  4b/ ?  /  1  C/R  x3dx  (37) 

“  TT  f\ 


COMPRESSIBILITY  AND  STALL 


The  expression  used  for  determining  profile  power  (6  vs  C.  )  does 

L(M) 

not  account  for  compressibi 1 i ty  above  drag  divergence  or  retreating  blade 
stall.  This  procedure,  using  the  mean  lift  coefficient,  does  not  compute 
the  variations  in  blade  angle  of  attack  or  local  lift  coefficients  along 
the  blade  or  around  the  disc.  To  account  for  the  effects  of  these  varia¬ 
tions  the  following  analysis  was  made. 


Compressibility  Above  Drag  Divergence 

It  is  assumed  that  the  unaccounted  for  compressibil i ty  effects  on  profile 
power  starts  at  1 /2 ( C^/a ) •  The  value  of  C^/o  at  the  start  of  drag 

divergence  was  termed  (C|_/o)^R,j.  The  following  method  is  used  to  calcu¬ 
late  (^l^^lsl* 

The  major  variables  which  contribute  toward  compressibil ity  effects  were 
assumed  to  be  u,  C^/a,  Cg/o,  and  eg.  The  condition  where 

bCqd/'->  =  .004  is  referred  to  as  the  "lower  stall  limit".  Analysis  of  test 
data  showed  that  for  any  given  Cg/o,  and  og  at  the  lower  stall 

limit  the  value  of  (1-m)/Cl/ oLSL  was  nominally  constant.  Utilizing  this 
relationship  yielded  a  means  of  estimating  the  onset  of  blade  stall. 

An  expression  relating  (1-m)/Cl/o)lSl,  Cq/o.  M(1>g())  and  oe  was  establish¬ 
ed  by  a  curve  fit  of  these  data.  (See  Figure  13  for  typical  cross  plot.) 
The  resulting  curve  fit  equation  relating  these  variables  is 

( ^l/c ) lsl  ~  ~M)/T(Cp/a,  ^]?90’  ee) 
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where 


f(CD/o,  M1j90,  ee)  =  [1750. (C0/a)2  -208.5  CD/a(3.6M2-3.71 (M+l ) ) 
+(((C D/o  +  .01)2/. 000133)-. 752)(M-.7)+38.7M-21M2-9. 9] 


(1.17+1.21750  ) 
e 


(M  =  M1 >go) 

Stating  that  the  increases  in  profile  power  due  to  compressibility  above 
drag  divergence  should  be  nominally  parabolic  and  a  function  of  the  operat¬ 
ing  and  critical  C./o.a  general  expression  for  the  power  increase  was  set 
up  as  L 

Cp  =  CT[f(CL/o-(CL/a)crit)]2/constant. 
comp  u 

The  function  equation,  f  (C, relating  the  operating  and 

critical  Cp/o  was  set  up  to  determine  the  difference  between  the  &  at  the 

operating  C.  and  the  increased  3  due  to  compressi Di 1 i ty.  These  func- 

L(M) 

tion  equations  were  evaluated  as  follows: 
for  u<.25 

f(CL/o-(CL/a)crit)=(l.-(.33+3.p)CL/o+.133M+.036)-(CL/o)crit<(39) 
for  c>.25 

f(CL/a-(CL/o)crit)=CL/o-(CL/o)crit-(1.08CL/o-.071)  (40) 


The  parabolic  constant  in  the  Cp  ^Q^p  equation  wac  evaluated  from  appropri 
ate  flight  test  data  as  0.01. 


Retreating  Slade  Stall 


When  the  operating  C^/o  increases  beyond  the  (C^/o)^9^(  retreating  blade 

5tali  occurs  and  a  further  power  demand  is  made  that  is  in  addition  to  the 
parabolic  (Cp^^p)  increase  discussed  above.  Test  data  'indicated  that  this 

increase  in  power  was  dependent  on  the  differential  oetween  the  operating 
C L/o  and  the  (CL/a)1 SL< 
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These  data  also  indicated  that  the  additional  power  increase  was  linear 
as  the  operating  C^/a  increased  beyond  the  (C^/g)^^^>.016->.020.  A  mean 

value  of  .018  was  used  to  establish  this  stall  divergence  point.  A  gener¬ 
alized  equation  for  computing  the  power  increase  due  to  the  stall  contri¬ 
bution  was  set  up  as 

aCP$  =  (C^/g-(C^/o )^$l )]°(C|yG“ (^l/ 0  - lsl^ 

When  C,/o  -  (C|_/°)LSL  <0.018  (41) 

f(CL/o-(CL/o)LSL)  =  3.0(CL/o-(CL/o)LSL) 
for  CL/o-(CL/u)LSL>  0.018  (42) 

f(CL/o-(CL/o)LSL)  =  833.33(CL/a-(CL/c)LSL)2 
+  30.65(CL/o-(CL/a)LSL)  +  .336 


Total  Horsepower  Correction 

The  total  horsepower  correction  due  to  compressibility  and  stall  was  then 
expressed  as 


Cp^  -  [f(CL/u-(CL/o)crit)V.01]oCT  +  [f(CL/o-(CL/u)LSL)]o(CL/a)LSL)a  (43) 


and  hp  =  Cp  PrrR2(nR)3/550 

s 


(44) 


The  coefficients  used  in  the  function  (f)  equations  were  derived  from 
flight  test  data  conducted  in  appropriate  flight  regimes. 


GROUND  EFFECT 


To  evaluate  rotorcraft  performance  when  operating  in  close  proximity  to 
the  ground,  a  means  of  accounting  for  the  ground  effect  had  to  be  estab¬ 
lished.  With  the  lack  of  substantiated  "classical"  methods  for  this 
analysis,  an  extensive  flight  test  program  was  conducted  with  an  $-66 
rotorcraft  at  Patuxent  River,  Maryland.  A  large  range  of  rotor  heights 
and  winds  were  investigated  during  these  tests.  These  data  were  nondimen- 
sionalized  with  respect  to  velocity  and  rotor  height.  It  was  assumed  that 
in  ground  effect  the  average  rotor  induced  velocity  is  reduced,  thus  re¬ 
ducing  rotor  induced  power  and  vertical  drag.  Tnerefore  these  data  were 
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incorporated  as  a  reduction  factor  of  the  required  induced  velocity,  the 
resulting  curves  from  this  analysis  are  shown  in  Figure  14.  The  symbol 
used  for  ground  effect.  A,  modifies  the  induced  power  term  to 


(6) 


The  value  of  A  was  set  equal  to  1  when  H/D  =  1,  as  power  differentials 
were  deemed  insignificant  above  this  height.  The  tests  conducted 
covered  lateral  airframe  protrusions  up  to  15t,  R  only.  The  validity  of 
these  curves  for  lateral  protrusions  in  excess  of  15,  R  has  not  been  sub¬ 
stantiated  . 

Using  these  curves  for  tandem  rotor  helicopters  the  roto'-  height  is 
assumed  to  be  the  average  height  of  the  two  rotors.  The  rotor  radius 
assumed  for  these  calculations  is  based  on  the  radius  required  to  obtain 
the  disc  area  for  a  single  rotor  which  is  equal  to  the  projected  area  of 
the  actual  system.  Reasonable  agreement  with  measured  flight  test  data 
has  been  obtained,  with  these  assumptions  for  tandem  rotors. 


Figure  14.  Ground  Effect  as  a  Function  of  Horizontal  Tip  Speed 
Ratio  and  Wheel  (SKID)  Height. 


DYNAMIC  FLIGHT 


The  previous  derivations  establish  a  means  of  calculating  the  steady-state 
power  requirements  for  a  given  rotor  system  and  operating  weight.  Con¬ 
versely  these  equations  can  also  be  used  to  calculate  the  thrust  developed 
by  a  given  rotor  system  absorbing  a  specified  power.* 

The  HDP-method  states  that  the  power  required  by  a  rotor  system  can  be 
computed  from  the  equation: 

hp  =  hp .  +  hp  +  hp  +  hp  +  hp  =  f'T) 
a  ri  ro  r  p  rc  K  s 

Due  to  complex  involvement  of  T  in  these  power  terms,  a  direct  solution 
for  (T)  is  impractical.  In  view  of  this,  an  iterative  procedure  can  be 
used  to  establish  the  value  of  thrust  which  satisfies  the  power  available. 
A  tolerance  of  0.5  horsepower  is  used  to  test  for  convergence. 

It  should  be  noted  that  the  tip  path  plane  angle,  aTpp,  the  horsepower 
available  hpa>  and  the  rotor  rpm  are  required  for  computing  the  thrust. 

CONTROL  INPUTS 


The  flight  path  a  rotorcraft  will  follow  is  a  function  of  the  orientation 
and  magnitude  of  the  thrust  vector  and  the  power  being  supplied  to  the 
rotor.  During  take  off  and  landing  maneuvers,  the  manner  in  which  the 
thrust  vector  is  manipulated,  through  cyclic  and  collective  input,  will 
define  the  flight  path  a  rotorcraft  will  follow.  In  order  to  compute  the 
space-time  relationship  of  a  rotorcraft  during  these  maneuvers,  a  means  of 
establishing  time  histories  of  the  tip  path  plane  angle  and  power  input 
must  be  determined.  This  analysis  assumes  that  the  thrust  is  perpendicu¬ 
lar  to  the  tip  path  plane  axis  and  the  tip  path  plane  axis  is  coincident 
to  the  control  axis  in  low  speed  flight. 


TIP  PATH  PLANE  ATTITUDE 


Although  ;he  pilot  has  complete  freedom  of  cyclic  input  control,  specified 
clearance  requirements  for  take  off  and  landing  maneuvers  demand  specific 
cyclic  input  rates  and  limits.  Flight  tests  conducted  for  confined  area 
operation  have  indicated  that  relatively  repetitive  inputs  for  a  specified 
procedure  are  required.  The  results  of  these  tests  are  shown  in  Figure  15. 
Analysis  of  these  data  indicated  that  the  time  histories  of  the  tip  path 
plane  attitude  could  be  normalized  into  two  equations.  These  equations 
can  be  expressed  for  the  tip  path  plane  transition  from  its  position  at 
the  start  of  an  event  to  the  maximum  angle  used  during  the  event.  The 


*  The  relationship  that  T  =  f (Gw)  is  only  applicable  to  steady  state 
flight  conditions. 
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second  portion  of  the  equation  covers  the  transition  from  the  maximum 
angle  used  to  the  angle  required  at  the  end  of  an  event.  These  two  transi¬ 
tions  can  be  defined  as  the  two  segments  of  an  event. 

An  event  is  defined  as  that  portion  of  the  maneuver  between  two  specific 
occurrences,  i.e.,  the  variation  in  tip  path  plane  attitude  between  lift 
off  and  the  engine  malfunction  point  ( EMP ) ,  would  be  the  event  between  the 
occurrences  of  Lift  Off  and  EMP. 

These  equations  were  expressed  for  the  two  segments  of  an  event  as 


t(K+l )-t 

aTpp(K+l)  =  ANGq  +  (ANGmax-ANGo)  sin  [t./2(- - - *)]  (45) 

RATE 

for  the  first  segment, and  for  the  second  segment 

V  -  V(K)  , 

aypp(K+1)  =  ujpp  sin  [n/2(-^-y - )]  +  tan  ’(D/T)  (46) 

where  subscript  K  is  the  time  index  counter 

0  is  the  condition  at  the  start  of  the  event 
Max  is  the  maximum  value  permitted 
Rate  is  the  maximum  rate  permitted 
X  is  the  desired  limit, 
z  is  the  phase  out  speed. 

Equations  45  and  46  yield  a  means  of  generating  a  time  history  of  tip  path 
plane  attitude  during  any  event.  It  should  be  noted  that  the  transfer  from 
equation  45  to  46  is  controlled  by: 

If  V(K)  <  zV^  use  equation  45 

If  t(K+l)  >  tRATE;  aTpp(K+l)  =  ANGMAX 


If  V(K)  >zVx  use  equation  46 
where  z  is  any  desired  fraction  of  Vx. 


Thus  any  tip 
fining,  ANG„ 


plane  time  history  can  be  generated  and  controlled  by  de~ 
V  ,  and  z  along  with  the  conditions  ANG_,  and  t_  at 


path 

"MAX’  tRATE’ 

the  start  of  any  event.  The  flight  tested  rates  and 
equation  45  and  the  results  superimposed  on  the  test 


o’  o 
limits  were  applied  to 
data  in  Figure  15. 


These  equations  of  tip  path  plane  attitude  yield  a  good  representative  time 
history  of  control  input  capability.  Since  this  normalization  of  tip  path 
plane  attitude  permits  an  infinite  number  of  combinations,  caution  should 
be  used  so  that  the  time  history  generated  are  within  reasonable  pilot  work 
load  limits  when  establishing  trends  and/or  takeoff  procedures. 
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POWER  INPUT 


Aside  from  the  tip  path  plane  attitude,  the  power  applied  to  the  rotor  and 
the  rate  at  which  this  power  is  applied  also  directly  affects  the  flight 
path  the  rotorcraft  will  follow: 

Power  (energy  per  unit  time)  can  De  supplied  to  the  rotor  in  three  ways: 
from  the  engines,  from  a  change  in  rotor  rpm ,  and  from  a  change  in  kinetic 
energy.  In  the  Helicopter  Dynamic  Performance  program  the  power  input 
from  .the  engine  is  generally  assumed  to  be  a  linear  function  with  time. 
Exceptions  to  this  linear  function  are  made  for  power  application  during 
landings  and  controlled  flight  path  takeoff  procedures. 

For  the  general  case  the  horsepower  available  at  any  time  point  can  be  ex¬ 
pressed  as 


$HP( K+l )  =  SHP(K)  +  HPMAX/(tRATE)  (47) 

up  to  the  specified  power  available  or  transmission  rating,  whichever  is 
lower. 

Power  application  rates  during  landings  are  determined  by  the  power  re¬ 
quired  to  maintain  the  desired  approach  rate  of  descent.  This  power  re¬ 
quired  is  determined  for  the  average  conditions  occurring  during  a  speci¬ 
fied  power  available  or  the  transmission  rating,  whichever  is  lower. 


ROTOR  ENERGY 


Rotorcraft  are  designed  to  prevent  wide  excursions  in  rotor  RPM.  This  is 
accomplished  by  coupling  the  collective  stick  with  the  fuel  control  or 
throttle  so  that  the  selected  RPM  is  maintained  with  collective  stick  in¬ 
puts.  When  the  maximum  engine  power  level  is  reached,  the  rotor  RPM  will 
change  with  further  increases  in  collective  settings.  Similar  to  the  Tip 
Path  Plane  variations,  the  pilot  has  full  control  of  these  changes.  There 
are  some  exceptions  to  this  which  will  be  discussed  later.  In  order  to 
evaluate  rotorcraft  capabilities  in  dynamic  flight,  the  effect  of  these  RPM 
variations  must  be  evaluated.  This  evaluation  considered  the  following 
relationships. 

The  power  that  is  released  or  absorbed  by  a  rotating  system  can  be  evalu¬ 
ated  by  considering  the  rotating  mass  of  the  system. 

The  stored  energy  in  any  rotating  system  can  be  expressed  from  the  general 
equation: 

ROTATING  ENERGY  =  [JM/2](n2) 

where  JM  is  the  mass  moment  of  inertia  of  the  rotating  system. 

If  a  change  in  the  rotational  speed  occurs,  the  change  in  energy  level  can 
be  expressed  as: 
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aROTATING  ENERGY  =  [d^K^2  -  «  2) 

By  specifying  the  time  interval  for  which  this  change  in  RPM  occurs  the 
•horsepower  be intf  released  or  absorbed  by  the  system  can  be  expressed  as 


HPfi-  =  [JM/1100](^2  -  U22)/At  (48) 

To  utilize  these  relationships  a  means  of  establishing  rotor  RPM  variations 
with  time  must  be  established.  Since  the  pilot  can  control  the  RPM  excur¬ 
sions  through  collective  adjustment,  the  basic  aerodynamic  drag-inertia  rela¬ 
tionships  are  inappropriate  for  this  analysis.  To  determine  typical  varia¬ 
tions  of  rotor  RPM  with  time,  flight  tests  were  conducted  under  saturated 
(max’num  power)  engine  conditions.  These  tests  were  conducted  for  specified 
collective  application  rates  and  RPM  excursion  limits.  The  results  of  these 
tests  are  shown  in  Figure  16.  These  test  results  indicate  that  for  the  spec¬ 
ified  conditions,  variations  in  rotor  RPM  with  time  are  reasonably  repeti¬ 
tive.  Analysis  of  these  data  further  indicated  that  the  change  in  RPM  over 
a  specified  time  period  can  be  generally  expressed  from  the  equation 


dNR  -  AliR  sin 


0/2( 


t( K+l )  -  t 


LRATE 


-)] 


where  dNR  is  the  RPM  change 

aNr  is  the  total  RPM  excursion 

tRATE  1s  the  totd^  excursion  time 

tQ  is  the  time  of  collective  application 

t(K+l)  is  the  time  at  any  interval  between  tQ  and  t^-^ 

The  results  of  this  equation,  computed  for  the  test  conditions,  were  super¬ 
imposed  on  the  test  data  in  Figure  16. 

On  this  basis  the  RPM  at  any  point  during  any  given  time  period  can  be 
expressed  as 


NR  '  \  +  ANR 
when  t(K+l)  -t 


sin  0/2  ( 


t(K+l)-t0 

tRATE 


)] 


-  ^ATE ;  NR  ’  V 


(49) 


where  ND  is  the  RPM  at  the  start  of  the  event. 

K0 

Further  flight  tests  were  conducted  to  establish  the  effects  of  -andom 
rates  of  collective  application  between  a  specified  RPM  excursion.  These 
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CHANC-E  IN  Nr 


VERTICAL  LIFT-OFF  FROM  WHEELS  LIGHT  CONDITION 


Figure  16.  Typical  Rotor  RPM  Change  Time  Histories. 


tests  were  also  conducted  under  saturated  engine  conditions.  The  results 
of  these  tests  are  shown  in  Figure  17.  The  generalized  equation,  calculat¬ 
ed  for  the  overall  RPM  excursion,  was  superimposed  on  these  test  data  for 
reference.  These  data  show  that  large  variations  in  the  rates  of  change 
in  RPM  occur  when  large  variations  in  collective  application  rates  are 
used.  To  further  evaluate  these  data,  ti.e  rate  of  change  in  RPM  with  time 
vr  plotted.  These  results  are  shown  in  Figure  18.  Analysis  of  these 
curves  showed  that  independent  of  the  manner  in  which  the  RPM  varied 
across  the  time  period,  the  total  energy  change  during  the  total  time 
period  is  the  same.  This  can  be  seen  from  the  equations 


t  dANp 

:'nr  -  0 


and  •  [J„/2](4  -q|) 


On  this  basis  it  was  stated  that  the  energy  transferred  by  the  rotor 
system  can  be  established  between  any  specified  end  points  during  an  event 
using  equation  49  to  establish  the  time  history  of  the  RPM  during  the 
event. 

It  should  be  noted  that  while  this  equation  might  not  yield  the  exact 
time  history  of  the  RPM  variation,  as  can  be  seen  in  Figure  17,  the  overall 
effects  on  the  flight  path  will  be  insignificant  as  the  total  energy  trans¬ 
ferred  is  accounted  for.  This  assumes  that  reasonable  collective  applica¬ 
tion  rates  are  used  and  the  period  of  RPM  change  consitutes  a  small  portion 
of  the  overall  flight  path  being  evaluated. 

There  are  some  areas  where  tne  pilot  does  not  have  full  control  of  the  RPM 
variations.  These  include  rapid  collective  application  during  take  off  and 
during  a  specified  collective  application  time  delay  at  the  point  of  en¬ 
gine  mal function. 

The  RPM  droop  or  decay  during  power  application  will  vary  with  different 
engine  models  or  control  unit  designs.  The  amount  and  rate  of  droop  is 
oest  established  from  flight  tests  conducted  with  the  particular  engine 
configuration  a. id  allowable  control  system  rates  being  used.  The  amount  of 
rotor  droop  during  a  specified  collective  delay  at  the  engine  malfunction 
point  is  a  function  of  the  rotor  inertia  and  collective  setting  at  the 
time  of  malfunction. 


The  collective  setting  required  by  the  rotor  system  is  a  function  of  the 
power  being  absorbed  by  the  rotor,  which  is  in  turn  related  to  the  thrust 
or  C-j.  being  produced.  A  relationship  between  collective  setting,  e^,  and 

rotor  thrust,  Cy,  for  any  given  flight  condition  can  be  established  from 


blade  element  theory  as 


(T+1.5h2)-~  +  1 . 5 >  (1-0.5/) 

_  0  Q _ 

( 1  -/+ 2 . 25/ ) 
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a  can  be  referenced  to  6^  as 

a  -  0^^  -<f  ( 51  ) 

Defining  p  as  ^  sina^  and  as  the  tip  path  plane  angle  relative  to  the 
airstream,  then  the  inflow  ratio  \  can  be  defined  as 


Vs  i  n  a  -  u 


Where  the  value  of  ;i  can  be  established  from  Wald's  equation  using  the  re¬ 
lationship 


u  =  u/u„ u  =  uu„ 
o  o  o 


substituting 

,\  = 


u  y.5CT 
B 


(52) 


Assuming  constant  induced  velocity  it  can  be  stated  that 

=?  =  tan  Vv)  (53) 

These  equations  are  then  used  to  calculate  the  change  in  rotor  RPM  for  a 
specified  change  in  collective  stick  position,  as  follows: 


For  a  specified  6^  and  flight  condition  the  rotor  tnrust  coefficient  (Cy) 
can  be  calculated  from  the  equation 


C 


T 


(l-u2+2.25u4)g-l .5>(l-0.5u2) 

(1+1. 5u2) 


(54) 


With  C,  known  the  power  required  for  the  flight  condition  can  be  calculat¬ 
ed  from  the  basic  energy  equations.  The  difference  in  the  power  required 
and  the  power  available  must  be  accounted  for  by  an  increase  in  power 
level  or  a  change  in  rotor  RPM  when  the  maximum  power  available  is  reached. 
Relating  this  power  differential  to  the  rotating  energy  equation  yields 


hpreq  "  hpav  Hpn  noo 


(;  !  2' -;2  '/'•  t 


Solving  for  the  RPM  at  point  2  yields 


c_ 

"c 


1100  HP.  ,t. 
- li  ~  3 


1/2 


(55) 
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The  term  (a)  in  equation  54  is  the  lift  curve  slope.  In  fixed  wing  and 
blade  element  theories  this  term  takes  on  values  of  5.73  with  variations 
that  depend  on  airfoil  section  and  element  Mach  numbers.  This  term  in 
rotorcraft  momentum  theory  corresponds  to  the  mean  lift  coefficient  (C.  ) 

L(M) 

slope  of  the  total  rotor  disc.  Values  of  (a)  can  be  established  by  solv¬ 
ing  equation  54  for  (a)  substituting  in  known  values  for  Cy  and  re¬ 
quired  for  the  initial  conditions  and  the  program  will  compute  the  lift 
curve  slope.  If  the  at  the  initial  point  is  not  inpu^the  program 

will  use  the  value  of  2.66  as  typical. 


DETERMINATION  OF  SPACE-TIME  RELATIONSHIPS 


The  previous  sections  established  a  means  of  determining  the  thrust  a 
rotor  is  capable  of  developing  for  any  given  power  input.  A  means  of 
establishing  rotor  RPM  levels,  tip  path  plane  orientation  and  horsepower 
levels  during  any  given  event  were  also  defined.  Considering  that  take¬ 
off  and  landing  performance  is  nothing  more  than  a  specified  series  of 
events,  then  the  flight  path  a  rotorcraft  will  follow  can  be  established 
by  specifying  the  occurrences  during  each  event  of  the  maneuver. 

This  analysis  establishes  the  flight  path  using  a  step  function  integral 
approach.  This  approacn  assumes  that  the  RP;-i ,  tip  path  plane,  horsepower, 
and  velocities  are  constant  during  any  given  time  period.  These  values 
are  based  on  the  computed  level  for  the  mid-point  of  the  time  interval 
used,  Witn  these  values  established,  for  the  time  interval  in  question, 
tne  tnrust  the  rotor  is  producing  can  be  determined. 

With  the  rotor  thrust  established,  the  unbalanced  forces  acting  on  the 
rotorcraft  during  the  time  interval  can  be  determined  from  the  equations 


TcoSaTpp 

1  KtL0NG 

sin  Hrpn  -D„ 

,nLAT  * 

(56) 

-  Tsin.i-j.pp 

LONG 

COS  a  -D 

,fklat  y 

(57) 

—  TeCS  Ol-jrpp 

‘^LONG 

COSaTPp  -(GW+Dw) 

ikklat  v 

(58) 

The  displacement  of  the  rotorcraft  can  then  be  determined  from  the 
equations 


a  =  F/M 

VNEW  =  V0LD  +  adt 
•:S  =  V0Lpdt  +  1/2  adt2 


(59) 

(60) 

(61) 
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S  =  r  ^ 

^TOTAL  o 

(62) 

These  basic  equations  of  motion  are  applied  for  the  unbalanced  forces  in 
the  X  (longitudinal),  Y  (lateral),  and  Z  (vertical)  directions  to  establish 
the  flight  path.  The  velocity  terms  in  these  equations  are  referenced  to 
ground  speed  using  the  relationship  -W^  y 

The  subscript  (  )  refers  to  the  wind  components  in  the  X,  Y 
tions  and  can  be  computed  from  the  equations: 

and  Z  direc- 

W(X)  =  Vwcos(DIRx-HDG)  cos(DIR^) 

(63) 

W^  =  Vwsin(DIRx-HDG)  cos  (DIRZ ) 

(64) 

W(z)  =  Vwsin(DIRx) 

(65) 

Where  DIR  is  the  wind  direction 
HDG  is  the  A/C  heading 
Vw  is  the  resultant  wind  velocity 

When  accelerations  are  imposed  on  an  object,  a  change  in  kinetic  energy 
occurs.  The  energy  expended  to  effect  this  change  can  be  expressed  in 
terms  of  horsepower  as: 

HPac(.  =  AKE/550  at  (66) 

Where  a KE  =  KE?  -KE] 

The  kinetic  energy  of  a  body  is  computed  from  the  equation 
KE  =  l/2MVa2 

Where  M  is  the  mass  of  the  body,  Va  is  the  absolute  scalar  velocity. 

This  acceleration  power  is  included  in  the  horsepower  available  term  (hpa) 

for  the  determination  of  rotor  thrust.  Since  the  accelerations  are  a  re¬ 
sult  of  the  time  step  calcul ation .the  following  process  is  used  for  evalu¬ 
ating  each  time  step.  Initially  the  acceleration  from  the  previous  time 
step  is  assumed.  For  the  first  time  step,  zero  acceleration  is  used. 

Based  on  this  assumption  the  unbalanced  forces  and  resulting  accelerations 
are  computed.  The  initial  and  resulting  accelerations  are  then  compared. 

2 

If  the  initial  and  resulting  accelerations  are  within  0.3  ft/sec 
(.0093  G's)  in  the  horizontal,  lateral  and  vertical  directions  the  time 
step  is  considered  to  be  converged.  If  these  tolerances  are  not  met  the 
initial  accelerations  are  set  equal  to  the  resulting  accelerations  and  the 
time  step  is  repeated  until  the  tolerances  are  met. 


CONTROLLED  FLIGHT  PATHS 


It  has  been  stated  previously  that  the  flight  path  a  rotorcraft  will  take 
is  entirely  dependent  on  tne  tip  path  plane  attitude,  the  rate  of  power  ap¬ 
plication,  and  the  maximum  power  used  during  the  maneuver.  If  the  flight 
path  a  rotorcraft  is  to  follow  is  controlled,  as  would  be  required  if  a 
height-velocity  envelope  existed,  then  the  rate  of  power  input  or  tip  path 
plane  attitude  must  be  controlled. 

Evaluations  of  height-velocity  diagrams.  Reference  7  resulted  in  a  means 
of  normalizing  the  boundaries  of  the  low  speed  height  velocity  envelope. 
This  normalized  diagram  is  shown  in  Figure  19.  Only  the  lower  portion  of 
this  curve  is  shown  herein. 


'  '  "  ■  GROUND  "J INC  '  - 


Figure  19.  Nominal  HV  Lower  Limb. 


’Using  these  data  the  actual  height  and  speed  at  any  point  along  the  curve 
can  be  determined  by  specifying  V  ,  h  and  y.  from  the  equations. 

A  A  i 

v  ■  <v/vx)vx 

and  Y  =  h1  (hx  -  y] )  +  y] 

Reference  7  specifies  the  coordinates  of  this  curve  in  tabular  form.  Analy¬ 
sis  of  this  data  showed  that  this  curve  could  be  reasonably  approximated  by 
the  equation 


Pegg,  R.  0. ,  AN  INVESTIGATION  OF  THE  HELICOPTER  HEIGHT-VELOCITY  DIAGRAM 
SHOWING  EFFECTS  OF  DENSITY  ALTITUDE  AND  GROSS  WEIGHT,  NASA  TN  D-m536 
May  1968. 
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A  comparison  of  the  tabulated  data  from  Reference  7  and  this  equation  is 
also  shown  in  Figure  19. 

This  equation  then  yields  a  means  of  defining  any  flight  path  by  specifying 
the  end  point  of  the  curve,  V  ,  h  and  the  initial  height  y, . 

In  order  to  establish  the  power  input  rate,  required  to  follow  this  flight 
path,  the  following  manipulations  of  the  h-|  equation  were  conducted. 

The  slope  of  the  curve  at  any  point  can  be  defined  by  taking  the  first  dif¬ 
ferential  of  h  with  respect  to  V,  thus 


h  =  y-y1  = 


dh/dv  = 


(1.1VX  -  V)2 


Multiplying  the  left  side  of  the  equation  by  dt/dt  yields 
[dh/dv] (dt/dt )  =  (dh/dt )/ (dv/dt )  =  Vy/aH 


therefore,  the  permitted  vertical  velocity  at  any  point  along  the  flight 
path  can  be  established  from 


V 


v 


[.nvx(hx-  yi) 
(i.ivx-v)2 


(68) 


Thus,  the  power  required  to  sustain  the  permitted  vertical  velocity  can  be 
established  for  the  time  interval  in  question.  It  is  assumed  that  the  tip 
path  plane  attitude,  ajpp,  is  specified  for  the  time  intervals.  It  should 

be  noted  that  if  the  computed  rate  of  power  application  is  greater  than 
the  maximum  allowable  rate,  the  maximum  allowable  rate  will  be  used,  also, 
if  the  power  requirements  exceed  the  maximum  allowable  power,  the  maximum 
allowable  power  will  be  used. 


Flight  paths  can  be  controlled  in  this  manner  (within  the  imposed  limits 
of  tip  path  plane  attitude,  power  application  rates  and  power  limits)  for 
any  predefined  flight  path.  The  flight  path  can  be  predefined  by  specify¬ 
ing  the  three  coordinate  points  y, ,  h  ,  and  V  . 

1  X  X 


5b 


Another  form  of  a  controlled  flight  path  would  be  to  specify  a  level  ac¬ 
celeration  at  maximum  power  up  to  a  given  speed  or  distance.  In  this  case 
the  power  is  allowed  to  vary  as  a  linear  function  up  to  the  maximum  avail¬ 
able,  or  transmission  limit,  and  the  tip  path  plane  attitude  required  to 
maintain  level  flight  is  computed.  The  tip  path  plane  attitude  is  com¬ 
puted  on  the  basis  that  if  level  flight  is  to  be  maintained  then  Tcos  aypp 

=  GW  +  D  v  or 


.  r  -1 /GW+D ,a 
aTPP  "  cos 


.where  T  =  f(SHP  &  a-rpp) 

This  requires  an  iterative  solution.  Caution  should  be  used  when  using 
this  type  of  controlled  path  to  insure  that  sufficient  ground  clearance  is 
maintained  as  high  tip  path  plane  attitude  requirements  could  result  from 
high  power  available  to  power  required  ratios. 
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AUXILIARY  POWER  DEVICES 


The  HOP  program  was  designed  to  accommodate  various  program  options  to 
provide  investigation  of  various  auxiliary  power  devices,  dropping  of  ex¬ 
ternal  loads  or  activating  drag  brakes.  The  various  available  devices  and 
the  equations  used  for  their  derivation  are  listed  below. 

1.  Auxiliary  Power  Units.  These  units  are  assumed  to  deliver  a  specified 
horsepower  to  the  main  rotor  gearbox  for  some  finite  time  period. 

This  power  is  added  to  the  horsepower  available  (if  any)  for  the 
determination  of  the  thrust  during  the  time  span  for  which  it  is 

appl ied. 

2.  Rotor  Tip  Rockets.  Tip  rockets  are  handled  by  converting  the  rocket 
thrust  to  horsepower  for  inclusion  into  the  power  available  term. 

This  rocket  thrust  is  converted  using  the  equation  hp  =  Tn  R  fi  b/550. 

dR  k 

The  HDP  program  assumes  that  if  tip  rockets  are  being  used,  there  is 
one  mounted"  on  each  rotct  blade.  Here  again  the  point  in  the  flight 
path  at  which  the  rockets  are  to  be  ignited  and  the  time  duration  for 
which  they  are  to  be  used  are  specified 

3.  Flywheel s .  The  use  of  flywheels  in  the  drive  train  is  handled  in  the 
same  manner  as  the  rotor  energy.  The  horsepower  contribution  to  the 
system  is  computed  using  equation  48, where  is  the  mass  moment  of 

inertia  of  the  flywheel  and  the  angular  velocity  is  computed  from  the 
relationship  P-MR  RPM  flywheel/RPM  rotor. 

If  flywheels  are  being  used  the  program  assumes  them  to  be  on  line 
during  the  entire  flight. 

4.  External  Propulsion  Units.  These  units  are  considered  as  JATO  or 
RATO  units  attached  to  the  airframe.  The  program  assumes  chat  these 
units  are  mounted  so  that  the  resultant  thrust  vector  passes  through 
the  airframe  center  of  gravity.  The  program  handles  these  thrust 
units  by  applying  the  horizontal,  lateral  and  vertical  components  of 
thrust  to  the  summation  for  determining  the  unbalanced  forces.  The 
basic  inputs  required  are  the  portion  of  the  flight  that  units  are  to 
be  applied,  the  thrust  per  unit,  the  number  of  units  being  used,  the 
mount  angle  of  the  unit,  and  the  time  duration  that  the  thrust  is 
available. 

5.  Dropping  of  External  Load.  Provisions  have  been  made  in  the  HDP 
Program  for  dropping  an  external  load,  if  applicable,  in  case  of 
emergency.  The  program  also  allows  for  a  change  in  airframe  drag  if 
this  option  is  used.  The  program  considers  this  delta  weight  change 
when  summing  forces  to  obtain  the  unbalanced  forces.  The  change  in 

1 F *  is  considered  in  the  evaluation  of  the  power-thrust  relationship. 

6.  Drag  Brakes.  To  minimize  the  program  code  the  use  of  drag  brakes  is 
handled  using  the  same  program  logic  block  as  used  for  dropping  an  ex¬ 
ternal  load. 
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In  this  case  the  load  dropped  would  weigh  0.0  pounds  and  the  change 

2 

in  * F '  would  be  +aF  ft  .  This  effectively  provides  evaluation  of  the 
drag  brakes.  In  both  cases,  items  5  and  6,  the  portion  of  the  flight 
during  which  the  load  is  to  be  dropped  or  the  drag  brakes  applied  can 
be  specified.  Once  this  option  has  been  turned  on  the  program  assumes 
that  it  applies  until  the  end  of  the  flight. 

The  auxiliary  power  devices  discussed  in  this  section  can  be  accessed  in 
the  HDP  program  through  input  locations  207-218.  Refer  to  Volume  II  for 
the  explicit  input  requirements  for  their  activation. 
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AIRFRAME  ATTITUDE 


The  treatment  of  airframe  pitch  and  roll  attitudes  in  the  HDP  program  is 
based  on  a  summation  of  moments  about  the  rotorcraft  center  of  gravity. 
Referring  to  Figure  20  it  can  be  seen  that  a  summation  of  moments  about  the 
CG  would  yield 

zM  =  _TTRDTR  +  HM  sln(SHAFT  TILT)  +  TMRDR  cos  otyppS i ns 


If  the  rotorcraft  is  in  a  steady-state  condition 

XM  =  o 


and  p  =  sin_1(-TTRDTR  +  sin(SHAFT  TILT)/(Tp.RDRC0SaTpp) 

and  AIRFRAME  PITCH  ANGLE  =  8  -Y 


Considering  the  rotorcraft  in  a  dynamic  state  with  a  time  variant  tip  path 
plane  attitude,  the  airframe  attitude  can  be  related  to  the  tip  path  plane 
attitude  by  considering  the  airframe  as  a  pendulum  with  a  pivot  point  at 
the  main  rotor  centroid. 


On  this  basis  the  unbalanced  moment  acting  on  the  airframe  can  be  expressed 
as 


PM  TTRDTR  '  (2  °R  TMR  +  'V  s1n^3B) 


(70) 


Where 


6  =  aTPPt  +  nt-l  "0o 


(71) 


Transferring  the  airframe  inertia  to  the  main  rotor  centroid 

I  =  I  +  w/gDD^ 
yyMR  yyCG  R 


(72) 


and  the  pitching  acceleration  would  be 

Prt  =  PM/I  =  l 
A  M  yyMR 


(73) 


The  airframe  pitch  velocity  during  the  time  interval  in  question  can  be 
expressed  as 


6t  *  Bt-i  +  6tat 


(74) 


59 


and 


(75) 


=  ^tAt  +  +  st-l 

It  should  be  noted  that  the  HDP  program  does  not  evaluate  airframe  aero¬ 
dynamics  be.yond  the  parasite  power.  Thus  airframe  lift  forces  and  pitching 
moments  are  not  included  in  the  abO'e  analysis.  To  account  for  these  items 
a  damping  term  was  applied  to  the  pitch  velocity  evaluation.  Thus  the 
pitch  velocity  equation  was  written 

Bt  =  +  ryt)PDA.iP  (74m) 


On  this  basis  PDAMP  (pitch  damper)  takes  on  the  definition  that  infinite 
damping  is  attained  when  PDAMP  =  0  and  no  damping  occurs  when  PDAMP  =  1. 
Thus  various  values  of  damping  can  be  evaluated  by  specifying  values  of 
PDAMP  between  0  and  1 . 

Tne  use  of  this  term  does  require  the  user  to  establish  the  value  of  PDAMP 
for  the  rotorcraft  being  evaluated  by  correlation  of  body  attitude  vs  time 
with  flight  test  measured  time  histories  of  dynamic  flight.  Once  this 
value  has  been  established,  trade-offs  can  be  conducted  to  evaluate  the 
effects  of  changes  in  static  stability  by  using  variations  of  PDAI1P  from 
the  base  value  and  ratio  the  PDAMP  values  with  the  ba^e  static  stability 
tc  find  the  new  static  stability. 

The  evaluation  of  airframe  roll  attitude  is  conducted  in  the  same  manner  as 
the  airframe  pitch  analysis.  When  dealing  with  the  roll  axis  the  airframe 
inertia  is  based  on  the  XX  or  longitudinal  airframe  axis.  The  basic  sketch 
showing  the  dimensions  and  forces  for  the  roll  analysis  is  shown  in  Figure 
21. 

Yaw  attitude  is  not  considered  in  the  HDP  Analysis. 

A  dynamic  body  attitude  analysis  for  tandem  rotors  has  not  been  conducted. 
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Force-Moment  Diagram  for  Establishing  Airframe  Pitch 


H-V  ENVELOPES 


The  HOP  program  has  incorporated  the  capability  of  computing  the  low-speed 
HEIGHT-VELOCITY  (H-V)  diagram.  When  the  program  is  operated  in  this  mode 
the  user  need  not  be  concerned  with  the  control  input  requirements  for 
program  operation.  The  program  is  designed  to  establish  the  best,  'optimum' 
flight  procedures  to  use  for  minimizing  the  avoid  area  of  the  H-V  enve¬ 
lope.  The  basic  input  requirements  for  this  program  mode  is  summarized  in 
section  6.3  of  the  User's  Manual  (Volume  II). 

The  construction  of  the  H-V  envelope  is  based  on  the  following  procedures. 
Reference  7  outlines  a  procedure  and  normalized  equations  for  estimating 
HEIGHT-VELOCITY  envelopes.  Using  this  as  a  basis.the  three  key  points  of 
the  H-V  curve  can  be  calculated. 

These  three  key  points  being  low  hover  height,  nose  point,  and  the  high 
hover  height.  Based  on  Reference  7,  the  low  hover  height  (h.  )  can  be  es¬ 
timated  using  the  equation: 

hl°  =  “  TV 00 "HP  A  d-2-24/0!-/0) 


To  account  for  partial  power  conditions  (multi-engine  rotorcraft)  the 
HP  A  term  was  modified  to  an  unsupported  power  loading  term.  Thus,  to 
account  for  any  partial  power  condition,  the  HP_A  terms  can  be  expressed 
as 


HP 


UNS 


hp.aO-hp^/hpjO 


and 


1100  HP 


UNS 


(1-2.24v/Ct/o) 


(76) 

(77) 


It  should  be  noted  that  if  the  horsepower  available  is  equal  to  or  greater 
than  that  required  to  hover  i.e.,  (HP^/HP^a  >1)  ,  then  no  HEIGHT-VELOCITY 

diagram  exists  and  further  evaluation  is  unnecessary. 

Using  this  estimated  h]g  point  as  a  starting  point,  the  program  cycles 
tnrough  a  double  iteration  scheme,  adjusting  the  h-  height  sc  as  to  satisfy 
a  user  specified  vertical  touchdown  velocity  and  c|^0 limit.  The  results  of 
this  iteration  yields  the  flight  procedures  required  and  the  maximum  low 
hover  height  from  which  these  imposed  conditions  can  be  met.  (For  further 
optimization  discussion  refer  in  text  to  the  Appendix.) 
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The  estimate  for  the  nose  point  critical  speed  is  outlined  in  Reference  7 
in  graphical  form.  This  figure  is  shown  as  a  function  of  velocity  for 
minimum  power  and  C^/o.  The  term  C^/o  is  in  turn  defined  in  terms  of  the 

minimum  power  speed  as 

CL/a  =  2/umin  Cy/ o  (78) 

Expressing  the  family  of  curves  given  in  Reference  7  jn  equation  form, 
yields 


VCR  =  2.809Vmin  +  5.618Cl/o  -169.776 

The  minimum  power  speed  (V^)  can  be  obtained  by  setting  the  first  dif¬ 
ferential  of  the  basic  power  equation  to  zero.  Expressing  the  basic 
speed-power  curve  as 


CT?  .  F  2  ,  06 


Cn  =  — 3-  + - 7  u  +3-  ( 1  +nMt' )  +  CTuy 

4  2MB  2  rr  8  1 


2 


Assuming  np  and  Cypy  (CLIMB  POWER)  small  at  minimum  power  speed  then 

Cn 


"T\2  1 


2  2 


3 Cq/ 3p  =  -1  /c(~r  ~  +  3/2  F/ n R  v 


B 


at  minimum  power  OC^/ap  =  0 
Solving  for  p  yields 


UM I N  ~  «VB>2  / 


To  better  approximate  the  minimum  power  speed  because  of  the  terms  dropped 
a  factor  of  0.95  was  used. thus 


UMIN  ~  0-95  ' 6F  ^  ^ 


(79) 


To  account  for  partial  power  conditions  the  critical  speed  (V 


CR; 


was  modi¬ 


fied  by  the  unsupported  power  loading.  For  this  condition  the  power  re¬ 
quired  to  hover  was  reduced  by  the  climb  power  differential  for  the  desired 
impact  speed,  or 
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(80) 


/  HP 

VCR  =  (2.809Vmin  +  5.618  ZJo  -169.776)  (,1 - 

(hr,  a 

Here  again  if  the  horsepower  available  is  equal  to  or  exceeds  that  required 
at  touchdown  no  HEIGHT-VELOCITY  diagram  exists. 

Based  on  the  results  given  in  Reference  7,  the  height  for  the  nose  point  is 
used  as  95  feet.  For  partial  power  conditions  available  test  data  indicate 
that  50  feet  is  more  appropriate  for  use  as  a  starting  estimate. 

The  high  hover  height  as  given  in  Reference  7  can  be  estimated  from  the 
equation 

n,  .  -  0. 18VrD2  +  199.0  (81 ) 

til  LK 


for  total  engine  malfunction. 

When  considering  partial  power  conditions  available  test  data  indicate 
that  the  relationship 


m 


GW_ 

nR* 


(.- 


HP 


AV 


HP,  A+ 


GW  V^ 


-)  * 


111.0 


(82) 


Yields  a  closer  approximation  to  actual  test  conditions. 

Evaluation  of  the  non  dimensional  shape  of  the  HEIGHT  VELOCITY  envelope  as 
given  in  Reference  7  indicates  that  the  following  relationships  can  be  used. 
For  the  lower  limb  (from  hj  to  nose  point)* 


\  =hi0  +  (O.n/O. i-vx/vCR)-.ii)(htR  -h|0) 

and  for  the  upper  limb  (h(1-  to  nose  point)* 

hx  =  hhi  ‘  0-(1-vx/vCr)  /Z)  (hhi  '  hCR^ 


The  equations  outlined  in  this  section  are  used  in  the  HDP  program  for  es¬ 
timating  the  starting  points  for  the  optimization  of  the  flight  procedures 
and  iterating  the  three  critical  points  to  attain  user  specified  ground 
contact  conditions.  The  HEIGHT-VELOCITY  curve  is  then  formed  by  using  the 


(83) 


(84) 


i 

t 
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PROGRAM  LOGIC  FLOW  CHARTS 

The  equations  contained  in  this  report  were  assembled  into  a  series  of 
subroutines  to  form  the  Helicopter  Dynamic  Performance  program.  A  brief 
summary  for  each  of  these  subroutines  is  given  in  Table  1.  The  subrou¬ 
tines  summarized  in  Tables  2  and  3  pertain  to  the  plotting  interface  rou¬ 
tines  provided  for  linking  with  the  TEKTRONIX  Plot  10  and  AG  II  plot  pack¬ 
age  and  the  Houston  Plot  Package  supplied  by  the  Government,  respectively. 

Program  logic  flow  charts  for  the  Engineering  analysis,  that  is,  all  routines 
contained  in  Table  1  have  been  provided.  Due  to  the  size  of  the  logic  flow 
charts,  these  charts  are  contained  in  a  separate  reference  compendium.  A 
copy  of  this  compendium,  in  microfisch  form,  may  be  obtained  from  ATL,  Aero¬ 
mechanics  technical  area  by  contacting  Mr.  T.  White,  phone  (804)  878-3874- 
-2062  or  autovon  927-3874/2062. 

While  computer  plots,  both  interactive  with  the  GT-40  scope,  TEKTRONIX 
4010,  4012  and  4014  and  Houston  plot  bed  capability  are  a  contractural 
requirement,  the  details  of  each  of  the  subroutines  are  beyond  the  scope 
of  this  contract.  The  basic  plot  package  consists  of  188  subroutines. 

It  will  suffice  to  say  that  the  interface  between  the  Engineering  Branch 
of  the  program  and  the  plot  package  interface  routines  provides  all 
parametric  plotting  capabilities  required  by  contract. 
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Table  1.  Concluded 
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5 

s 

O 

O 

GH 

UJ  u. 

UJ  Q 

UJ 

0 

3 

GB 

GE 

GC 

GF 

a  z 

ca 

-j 

c 

-  . 

• 

o 

1“  £: 

p—  • 

x 

£ 

tO  — •  3 

*£  Z  CQ 

ja— -< 

4-> 

HZO 

-J 

to 

o 

0 

O 

Q. 

♦— 

H- 

CQ  >-  X 

X 

3  . 

—I  'X  CO 

a 

O 

h- 

! 

—  x  < 

UJ 

UJ 

— J 

< 

o 

O 

Q. 

£  CQ  CQ 

CQ 

*“  H-  — 1  O 

•  ce  £  < 

to 

i/o  —  o  £ 

>  to 

>*  to 

Z 

J  U.  X  _J 

LU  Z 

O  uj 

O  UJ 

b-  X  to  X 

CD 

oc  *-«  *  *• 

re  o 

QC  UJ 

Q 

O'  r-  r- 

Q-  z 

CL  Z 

X  •  x 

UJ 

UJ 

3^ 

►—  O  Z 

>- 

UJ  QC 

UJ  QC 

►  ruxo 
zSau  j 

X 

to 

1 

0) 

*  -J  O  UJ 

£ 

Z  UJ 

Z  UJ 

O  _» 

£ 

£ 

1 

to  cl  £  to 

— •  U. 

•— «  u_ 

QC  CQ 

— •  u.  X  UJ  CL 

O 

1 

|  C  I 

UJ  Z  Z  H- 

-J 

b-  UJ 

b-  UJ 

<  £ 

-J  £  b-  Z  t_J 

_i 

O 

ao< 

to 

3  QC 

3  QC 

CO  >- 

QC  X  X  O 

to 

OX 

O 

O 

to 

£  -  .  -* 

p 

X 

UJ  <— 

QC  >- 

QC  >- 

r-  O 

U. 

Lt_ 

to  * 

—  X 

b-  X 

3 

•  a  uj  < 

X 

>-  tO 

» —  ■ — 

*•  £  uj  tO  C 

X 

UJ  UJ 

£  £- 

«—  UJ 

►  -  QC  Q.  CL  -J 

£ 

5?» 

o  >■  <t  a 

> — 1 

£  Z 

_ 1  U.  X  UJ  CL 

_J  ►— «  I  UJ 

Z  3  < 

3  <£ 

Q.  £  b-  b-  tO 

O  CO 

-O 

DCOQ.Z 

to 

£  owi 

O  CO 

_)  _1 

UJ  X  X  to  O 

to 

(_>  — . 

rx 

oo 

<U 

u 

! — 

5 

o 

0 

CL 

ro 

b~ 

H- 

X 

O 

UJ 

UJ 

5- 

-J 

-J 

CL 

<u 

CQ 

o 

<c 

•a; 

+-> 

o 

o 

C 

h- 

_J 

_j 

z 

to 

Q. 

Cl 

0 

CL 

< 

c 

X 

X 

O 

CL 

h— 

►— 

0 

z 

UJ 

_j 

o 

0 

CO 

z 

0 

c 

O. 

o 

z 

UJ 

3C 

o 

0 

i- 

* 

O 

0 

0 

* 

UJ 

UJ 

<_> 

o 

O  oc 

to 

X 

o 

~J 

—1 

QJ 

Q-  J — 

CL  H* 

z 

cl  o 

CL  t_> 

1 - 

UJ 

UJ  3  <t 

S3 

UJ 

0 

3 

Cl  tO  QC 

CL 

-J 

t 

Si 

o 

>-  b~ 

>- 

8 

<L> 

z 

h-  t~  Z 

• —  Z 

b- 

z 

ce 

0 

—1  o  o 

O  O 

.J 

O 

cn 

j—  ^ 

3 

Cl. 

a.  z  o 

z  0 

CL 

V- 

3 

*-* 

h- 

ns 

u 

rd 

a. 

►— 

z 

4-> 

UJ  »— 

o 

£  tO 

3  — 

r*j 

f'J 

1 

rsi 

1  i  ■ 

O  -J 

*> 

l 

i 

1 

» 

_ 

l-vl 

Q- 

< 

>- 

>- 

>- 

X 

* 

• 

« 

X 

X 

X 

X 

S- 

UJ 

a) 

o 

UJ 

to 

Z 

b-  to 

c 

< 

o 

O  UJ 

a» 

—J  z 

a_  — 

O 

CZJ 

< 

Cxi 

_j 

O 

5 

o 

x  (j 

0 

to 

>> 

— J 

O 

UJ  z 

z 

to 

o 

»— * 

z 

UJ 

to 

5 

UJ 

fsl 

g  s 

►- 

b- 

U. 

UJ 

O 

3 

i- 

O 

QC 

u.  o 

O 

o 

U- 

QC 

to 

to 

-J 

cs  a. 

H- 

b- 

CL 

O 

UJ 

z 

O 

O 

QC 

to 

O 

§ 

— •  to 

-J 

-J 

g 

O 

< 

O 

CL 

CL 

to 

00 

Cl. 

*— 

oc 

3  3^ 

U. 

o  o 

u. 

h- 

o 

to  Q.  UJ 

h- 

< 

a. 

o  z 

3 

< 

C 

H" 

UJ 

UJ 

Z  O 

>-  QC  _J 

Z 

to 

UJ  UJ 

o 

^  — 

—J  3  < 

U> 

a: 

a 

< 

C\J 

3  O  > 

< 

UJ 

<  3 

QC 

> 

X 

to  o 

QC 

►—  _J 

u.  z  h- 

a 

1 

3 

O  b- 

<u 

<  o 

UJ  05 

0  0  < 

cr 

Z  3 

oc 

r-J 

to 

CQ 

X 

UJ 

— *  CL 

3 

»-•  to 

3  UJ 

UJ 

to 

H 

XI 

UJ 

a. 

tO  3 

c  o 

to  •— < 

b~ 

3  3  Z 

0 

3 

3 

3 

h-  O 

to 

n3 

— «  O c 

3  to 

O 

to 

»—  »— 

o 

3  < 

to  to  0 

> 

to 

to 

5  0 

h— 

— 1  o 

3S 

b- 

H 

z  o 

mJ  Z 

UJ  UJ  Lu 

Ul 

UJ 

UJ 

UJ  «J 

X 

•—  o 

<  < 

<  UJ 

to  to  O 

Q. 

to 

to 

to 

CL  CL 

UJ 

to  ►— 

3 

Z 

_J 

>—  -J 

QC 

-J 

a. 

o 

O 

< 

Q. 

h- 

_J 

-J 

CO 

o 

QC  QC 

£  z 

> 

— 1 

QC 

»-* 

3  < 

3 

-J 

O'  O 

— *  0 

Z 

CL 

£ 

o  z 

DC 

H- 

3 

UJ  UJ 

CQ 

3 

UJ 

1— 

74 


(1)  TEKTRONIX  PLOT  10  ROUTINES  TO  BE  SUPPLIED  BY  THE  GOVERNMENT . 


Table  2.  Concluded 
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(1)  Tektronix  Plot  10  Routines  to  be  Supplied  by  the  Government. 


Sikorsky  General  Plot  Package-  Houston  Plotter  Interface-Subroutine  Description 
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(1)  TEKTRONIX  PLOT  10  ROUTINES  TO  BE  SUPPLIED  BY  GOVERNMENT 

(2)  HOUSTON  PLOT  ROUTINES  TO  BE  SUPPLIED  BY  GOVERNMENT  AND  MODIFIED  BY  SIKORSKY  AIRCRAFT 
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(2)  HOUSTON  PLOT  ROUTINES  TO  BE  SUPPLIED  BY  GOVERNMENT  AND  MODIFIED  BY  SIKORSKY  AIRCRAFT. 
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APPENDIX  A 


A  METHOD  FOR  EVALUATING  AUTOROTATION 

To  the  engineer  and  designer  a  "good"  autorotation  will  always  result  in  a 
"safe"  power  off  landing,  and  efforts  are  directed  toward  obtaining  the 
highest  possible  hover  height  for  the  lower  limb  of  the  height-velocity 
(H-V)  curve,  the  lowest  possible  hover  height  for  the  upper  limb,  and  the 
mini  mum  possible  critical  speed  for  the  "nose"  of  the  curve.  Practical  fac¬ 
tors,  such  as  the  ability  to  accurately  measure  airspeed,  or  concern  for 
control  input  and  field  size  requirements  are  generally  not  considered. 

Pilots'  comments  and  "feelings"  on  autorotation  capabilities  can  be  re¬ 
lated  to  rotor  and  rotorcraft  parameters  which  should  be  considered  in 
the  design.  A  selection  of  some  of  these  comments  and  the  related  tan¬ 
gible  parameters  are  listed  below. 

Pilot  Comment:  “Engine  malfunction  can  occur  at  any  time.  I  don't 

want  restrictions  placeo  on  my  flight  envelope." 

Tangible  Parameter:  Design  consideration  for  minn  the  avoid  area 

of  the  hei ght-veloci ty  diagram. 

Pilot  Comment:  "After  engine  malfunction,  there's  not  much  time 

until  I'm  on  the  ground." 

Tangible  Parameter:  Provide  low,  controllable  rates  of  descent.  Design 

considerations  for  disc  loading,  rotor  RPM  and  air¬ 
foil  type. 

Pilot  Comment:  "It  wouldn't  be  bad  if  I  were  over  Kansas  but  all  I 

^  can  find  is  a  small  rough  clearing." 

Tangible  Parameter:  Design  for  reasonable  roll  on  velocities,  landing  ge<  r 

design,  tai •  wheel  or  structurally  and  psychologically 
sate  tail  r,  '.or  protection. 

Pilot  Com  .r  "When  1  get  into  a  iiiyh  cyclic  flare  I  loose  all 

ground  re  arence. " 

Tangible  Parameter:  Cockpit  visibility,  minimize  attitude  change  require- 

pi  rv  p4-/* 

mu  i '  vj  • 

While  these  are  not  all  of  the  pilot  reactions  tc  emergency  autorotati on, 
they  are  pertinent  and  provide  a  starting  point  for  an  approach-  to  a 
meaningful  analysis  of  autorotati ve  characteristics 
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In  order  to  evaluate  the  effect  of  design  or  parametric  changes  on  the 
autorotative  capabilities,  a  means  of  extracting  the  maximum  capabilities 
for  a  given  parametric  set  must  be  established.  The  Helicopter  Dynamic 
Performance  (HOP)  program  provides  a  means  of  evaluating  rotorcraft  take¬ 
off  and  landing  performance  with  power  on,  and  with  partial  or  total 
engine  malfunction.  The  transition  phase  between  power  on  and  power  off, 
autorotative  entry,  is  also  included  in  the  HDP  program.  Using  this  pro¬ 
gram  as  a  base,  several  subroutines  were  written  and  inter-faced  with  the 
base  HDP  program.  These  subroutines  were  designed  to  generate  the  program 
inputs  required  to  minimize  the  height-velocity  (H-V)  envelope  for  a  given 
configuration  description  and  specified  end  point  constraints.  Based  on 
pilot  comments,  the  most  pertinent  end  point  constraints  were  selected  as: 


1.  vertical  velocity  at  touchdown 

2.  blade  loading,  Cj/c ,  at  touchdown 

3.  body  attitude  at  touchdown 

4.  horizontal  roll-on  speed. 

Based  on  these  end  point  constraints,  subroutines  were  written  to  deter¬ 
mine  the  point  in  space  where  the  rotorcraft  would  have  to  be  and  the 
flight  procedures  to  be  used  to  best  satisfy  the  imposed  end  point  limits. 
Three  routines  were  written  for  this  purpose.  One  for  the  low  hover 
height?  one  for  a  forward  speed  approach, and  the  third  for  the  high  hover 
height. 


Low  Hover  Height 


The  low  hover  height  routine  establishes  the  maximum  hover  height  from 
which,  if  a  total  engine  malfunction  occurred,  tne  rotor  craft  could  land 
without  exceeding  an  imposed  vertical  touchdown  speed  or  the  imposed  blade 
loading  (Cx/n).  This  routine  assumes  a  pure  vertical  descent  so  that  the 
roll  on  speed  and  the  body  attitude  criterio  are  not  pertinent.  Proaram 
output  for  this  condition  provides  the  maximum  low  hover  height  and  the 
collective  control  application  procedures  required  to  effect  the  landing. 
All  time  history  data  provided  by  the  main  program  is  also  available. 

Minimum  Speed 

The  second  routine  deals  with  a  total  engine  malfunction  occurring  during 
forward  flight.  The  purpose  of  the  routine  is  to  establish  the  minimum 
forward  speed  and  altitude  from  which  a  landing  can  be  effected  witnin 
the  imposed  end  point  constraints  The  flight  procedures  required  to 
effect  this  landing  are  also  a  program  output  requirement. 


Initially  it  would  appear  that  the  only  requirement  is  to  adjust  the  con¬ 
trols  from  the  point  of  engine  malfunction  to  touchdown  at  the  specified 
vertical  impact  speed  while  bringing  the  C T/i  and  body  attitude  to  the 
stated  limit  and  bringing  the  horizontal  roll  on  speed  to  the  lower 
specified  limit.  A  more  detailed  review  of  the  stated  end  point  con- 
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straints  shows  that  certain  combinations  of  innut  criteria  could  result 
in  conflicting  conditions.  For  example,  if  a  10-foot-per-second  vertical 
impact  speed  is  requested  and  a  maximum  Cj/o  of  0.20  is  used,  analysis 
for  some  disc  loading  conditions  could  show  that  the  rotorcraft  would 
have  a  brief  positive  rate  of  climb.  Another  conflicting  situation  could 
occur  if  a  low  body  attitude  at  touchdown  was  requested  along  with  a  low 
horizontal  roll  on  speed.  In  this  case  the  body  attitude  restriction 
could  preclude  the  deacceleration  required  to  reach  the  imposed  roll¬ 
on  speed  requirement. 

In  view  of  these  potential  conflicts,  the  end  point  constraint  criteria 
were  set  up  on  a  priority  basis  with  the  following  interpretations. 

1.  Vertical  touchdown  criteria. 

This  parameter  has  first  priority  and  must  be  met 
within  +0.5  foot  per  second  of  the  requested  value. 

2.  Blade  loading,  C./o. 

This  input  has  second  priority  and  is  interpreted  as  the  maxi¬ 
mum  allowable  value.  The  program  will  accept  a  lower  value 
if  the  specified  limit  is  not  consistant  with  the  other  con¬ 
straints  . 

3.  Touchdown  body  attitude. 

This  input  retains  a  high  priority  status  provided  the  tip  path 
plane  attitude  required  to  attain  the  permitted  touchdown  angle 
does  not  cause  the  rotorcraft  to  accelerate  during  the  final 
phase  of  the  maneuver.  If  acceleration  is  required  to  meet  the 
criteria,  the  tip  path  plane  angle  is  set  to  -1  degree  and  tne 
touchdown  body  attitude  limits  are  ignored.  The  program  will' 
print  out  a  diagnostic  message  if  this  condition  exists. 

4.  Horizontal  roll-on  speed. 

This  input  limitation  carries  the  lowest  priority  of  the  four. 
Since  the  main  purpose  of  these  subroutines  is  to  minimize  the 
height-velocity  envelope  and  the  other  three  criteria  carry  a 
higher  pilot  concern,  the  program  will  establish  the  flight, 
procedures  which  maximize  the  other  criteria  while  minimizing 
the  horizontal  approach  speed  and  roll-on  speed  independent 
of  the  roll-on  speed  requirement.  It  should  be  noted  that  the 
program  will  work  toward  getting  the  rotorcraft  within  the 
specified  criteria  so  that  the  inputs  for  these  values  should 
not  be  ignored.  The  program  will  print  out  a  message  if  the 
input  requirements  cannot  be  met. 
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High  Hover  Height 


A  third  routine  has  been  set  up  for  determining  the  high  hover  height 
of  a  height-velocity  envelope.  This  high  hover  height  can  be  defined 
as  the  minimum  height  at  which  a  rotorcraft  can  hover  and  effect  a  safe 
landing  after  total  engine  malfunction.  It  differs  from  the  low  hover 
height  in  that  acceleration  to  some  forward  speed  is  required  to  safely 
effect  the  landing.  The  same  end  point  constraint  criteria  and  priority 
order  and  definition  are  used  for  this  routine  as  for  the  minimum  speed 
point. 

Other  Constraints 


Aside  from  the  end  point  constraints  as  discussed  above,  the  program  user 
can  also  specify  the  maximum  range  of  cyclic  (tip  path  plane),  collective 
(THETA  75),  and  minimum  and  maximum  rotor  RPtl  within  which  the  program  can 
operate.  These  inputs  are  interpreted  as  the  maximum  limits  and  the  pro¬ 
gram  will  work  within  the  range  imposed,  but  not  necessarily  use  them 
if  a  lower  value  produces  a  "better"  flight  path.  The  "better"  flight 
path  is  defined  as  that  path  which  best  satisfies  the  end  point  con¬ 
straints. 

Sesi gn  Cri teri a  E val uation 


These  height-velocity  routines,  coupled  with  the  HOP  program,  provide  a 
means  of  evaluating  various  rotorcraft  design  criteria.  The  premise  for  this 
is  based  on  the  ability  of  these  routines  to  compute  the  minimum  boundaries 
of  the  H-V  envelope  for  a  given  rotorcraft  configuration,  control  limits, 
and  end  point  constraints.  Then,  by  the  perturbation  of  any  one  pertinent 
parameter,  the  resultant  change  in  the  H-V  envelope  is  a  measure  of  how  the 
safety  of  the  base  rotorcraft  is  improved  or  compromised.  A  systematic  para¬ 
metric  evaluation  can  then  highlight  the  design  parameters  which  hav»  signifi¬ 
cant  effects  on  autorotative  capability. 

Note  that  these  routines  generate  the  control  input  requirements  to  maxi¬ 
mize  rotorcraft  capabilities  without  regard  to  pilot  work  load.  The  H-V 
envelope  predicted  could  therefore  be  beyond  realistic  pi lot-rotorcraft 
capability.  See  Section  4,  Volume  II. 
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PARAMETRIC  EVALUATION  OF  AUTOROTATIVE  CAPABILITIES 

The  HOP  program  was  used  to  evaluate  two  helicopters  in  autorotative 
flight.  The  purpose  of  this  study  was  to  determine  what  effect,  if  any, 
various  design  parameters  have  on  the  autorotative  capability  of  these 
rotorcraft.  The  rotorcraft  used  for  this  study  were  specified  by  contract 
as  the  Hughes  0H-6A  and  the  Bell  AH-16.  The  parameters  to  be  varied 
are  listed  below: 

1.  Disc  Loading 

2.  Altitude 

3.  Center- of-Gravi  ty  Location 

4.  Rotor  Inertia 

5.  Maximum  Flare  Attitude 

6.  Rotor  Profile  Drag 

7.  Pilot  Control  Initiation  Delay  Time 

8.  Initial  Rotor  RPM 

9.  Aircraft  Parasite  Drag 

10.  Airframe  Static  Stability 

)i.  Vertical  Touchdown  Velocity 

12.  Touchdown  Body  Attitude 

13.  Maximum  Cy/c 

Three  additional  parameters  were  deemed  as  not  appropriate  for  the 
following  reasons: 

1.  Airspeed  -  The  airspeed  required  at  the  nose  point  of  the  H-V 
envelope  is  a  result  of  the  calculations  for  minimizing  the 
H-V  envelope  and  cannot  be  a  program  input. 

2.  Initial  steady- state  rate  of  climb  and  descent.  The  use  of 
approach  conditions  other  than  hover  and  steady-state  level 
flight  for  H-V  determination  is  not  pertinent. 

3.  Horizontal  touchdown  speed.  This  is  a  program  input.  Due  to 
touchdown  criterion  priorities,  the  use  of  this  criterion  for 
explicit  variations  in  capabilities  might  not  be  meaningful. 
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The  evaluation  of  each  rotorcraft  was  conducted  by  first  establishing  a 
baseline  condition  for  the  given  helicopter.  Aside  from  the  basic  air¬ 
craft  configuration  data,  baseline  conditions  of  sea  level  standard 
atmosphere  and  neutral  CG  were  used. 

Initial  calculations  were  conducted  to  establish  the  H-V  envelope  for  each 
rotorcraft  at  a  disc  loading  of  4.  The  results  of  these  calculations, 
shown  in  Figure  A-l(a),  show  that  significant  differences  exist  for  these 
aircraft,  in  autorotative  capability,  at  these  baseline  conditions. 

To  obtain  a  more  meaningful  evaluation  of  parametric  variations,  the 
perturbations  were  conducted  from  a  baseline  H-V  curve  which  would  pro¬ 
vide  a  common  H-V  speed  range  for  both  aircraft.  While  an  exact  common¬ 
ality  could  not  be  found,  a  reasonable  common  base  was  determined  with 
the  0H-6A  at  a  disc  loading  of  3.5  and  the  AH-1G  at  a  disc  loading  of  5. 

A  comparison  of  the  baseline  H-V  curves  is  shown  in  Figure  A-l(b)  para¬ 
metric  listing  of  baseline  conditions  is  in  Table  A-l  for  the  0H-6A 
and  Table  A-2  for  the  AH-1G.  The  parametric  variations  used  are  shown 
in  the  diagonal  or  staggered  boxes  in  Tables  A-l  and  A-2.  The  effect 
of  each  parametric  change  was  then  established  by  computing  the  minimum 
H-V  envelope  for  the  imposed  conditions  and  comparing  the  differences 
from  the  baseline  curve.  The  differences  between  the  H-V  curves  can  be  de¬ 
fined  as  the  change  in  low  hover  height,  approach  speed  at  the  "nose" 
point,  high  hover  height,  and  the  horizontal  roll-on  speed  from  the  "nose" 
point  and  hsgh  hover  height.  The  resulting  H-V  curve  and  variations  from 
the  baseline  are  tabularized  on  the  lower  portion  of  Tables  A-l  and  A-2. 
These  variations  from  the  baseline  are  also  shown  in  bar  graph  form  in 
Figure  A-2  for  the  0H-6A  and  in  Figure  A- 3  for  the  AH-1G.  (Sea  Level 
standard  conditions  apply  except  as  otherwise  noted.) 

A  review  of  these  figures  shows  that  while  the  sensitivity  of  each  para¬ 
metric  variation  differs  between  aircraft,  the  trend  of  increasing  or 
decreasing  the  H-V  envelope  is  nominally  the  same.  A  generalized  over¬ 
view  of  the  effect  of  each  parametric  change  is  presented  in  Table  A -3. 


DISCUSSION  OF  PARAMETRIC  TRADE-OFFS 


An  analysis  has  been  conducted  to  establish  the  effect  of  basic  design  para 
metric  changes  on  the  autorotative  capabilities  of  the  representative 
rotorcraft.  While  this  study  points  out  the  trends  of  each  parametric 
change,  it  also  points  out  that  the  magnitude  of  the  change  was  dependent 
on  the  aircraft  in  question.  Based  on  this,  it  cannot  be  said  that  a  cer¬ 
tain  parametric  change  is  going  to  have  a  given  effect  on  all  rotorcraft. 

An  example  of  this  C3 n  be  seen  with  the  effects  of  increasing  rotor  in¬ 
ertia.  Analysis  of  the  OH-6A  shows  that  a  50:1  increase  in  the  rotor 
inertia  decreases  the  overall  H-V  envelope  and  permits  a  lower  roll-on 
speed  from  both  the  "nose"  point  and  high  hover  height  condition.  On 
the  other  hand,  the  same  percent  increase  in  rotor  inertia  for  the  AH-1G 
raises  the  low  hover  height,  shows  no  change  in  the  "nose"  point  speed  and 
requires  a  slight  increase  in  the  high  hover  height.  The  resulting  H-V 
envelopes  for  this  evaluation  is  shown  in  Figure  A-4.  Review  of  the  roll¬ 
on  speeds,  Table  A -2  shows  a  slight  decrease  from  the  "nose"  point  and  an 
increase  from  the  high  hover  height. 

From  these  results  it  appears  that  the  sensitivity  of  the  rotor  inertia 
parameter  is  dependent  on  the  baseline  rotorcraft.  While  this  is  partial¬ 
ly  the  case,  one  must  keep  in  mind  the  priority  of  the  end  point  con¬ 
straints  imposed  on  the  analysis.  The  analysis  stated  that  the  ver¬ 
tical  impact  speed-  and  touchdown  body  attitude  have  top  priority  and 
blade  loading,  Cr/a,  arid  horizontal  roll-on  speed  have  second  priority. 
Based  on  this,  a'review  of  the  maximum  CT/o  levels  used  in  affecting 
the  landing  would  provide  further  insight  as  to  potential  effects.  A 
comparison  of  the  maximum  Cy/o  levels  required  to  effect  the  landing 
with  the  base  aircraft  and  rotor  inertia  changes  is  shown  below. 
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It  can  be  seen  that  for  the  01I-6A  the  maximum  C y/o  reached  is  nominally 
the  same  for  both  the  baseline  and  increased  inertia  case.  The  resulting 
H-V  envelopes  show  significant  improvement  with  the  increased  inertia. 

For  the  Ah'-IG  the  low  hover  height  shows  the  same  maximum  Cy/a  was  used 
and  an  increase  in  the  low  hover  height,  comparable  with  the  0H-6A,  was 
also  shown.  On  the  other  hand,  while  the  "nose"  point  and  high  hover 
height  maximum  C_/,:  was  near  maximum  for  the  baseline  AH-1G  (a  maxi¬ 
mum  value  of  0.1!>  was  imposed) ,  the  Cy/o  level,  used  with  increased 
inertia  had  to  be  reduced  to  meet  the  end  point  constraints.  From  this, 
it  can  oe  seen  that  if  the  AH-1G  pulled  the  same  Cy/;-.  with  the  i  creased 
rotor  inertia,  the  vertical  impact  speed  and/or  the  horizontal  roll -on 
speed  could  be  reduced.  It  should  be  noted  that  whi’o  this  roll-on 
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speed  reduction  appears  possible,  the  touchdown  body  attitude  limit 
might  be  the  constraining  factor  which  precludes  the  use  of  the  tip 
path  plane  angle  required  for  effecting  the  speed  reduction. 

Conversely,  it  can  be  stated  that  by  maintaining  the  baseline  H-V  enve¬ 
lope,  increasing  the  rotor  inertia  will  increase  the  level  of  safety 
for  specified  end  point  constraints. 

This  leads  to  the  conclusion  that  one  cannot  always  point  to  a  singular 
parameter  and  state  "this  is  the  one  that  will  make  things  better." 

As  in  the  case  discussed  here,  a  50G  increase  in  rotor  inertia  for  the 

0H-6A  shows  a  significant  reduction  in  the  H-V  envelope  while  using 
initially  the  same  blade  loading,  Cj/a  potential.  On  the  other  hand 
the  same  inertia  increases  for  the  AH-1G  shows  no  significant  change 
in  the  H-V  envelope  (except  the  low  hover  height)  primarily  because  the 
permitted  body  attitude  at  touchdown  does  not  allow  the  rotor  to  use 
the  available  blade  loading  potential. 

To  obtain  a  more  meaningful  comparison  of  parametric  effects  of  autorota¬ 
tion,  a  rating  factor,  based  on  the  ratio  of  blade  loading  potential  and 

blade  loading  used  for  attaining  a  given  H-V  envelops  is  introduced,  i.e., 

Cy/j  allowed/  C j/a  used.  This  rating  factor  has  a  baseline  value  of  1.0. 
The  efft-Ct.  of  any  given  parametric  change  is  basically  reflected  in  the 
change  in  the  H-V  envelope.  In  conjunction  with  this  H-V  change,  varia¬ 
tions  of  the  rating  factor  from  1.0  indicates  that  this  parametric  change 
could  be  more  effective  provided  it  was  changed  in  conjunction  with  an¬ 
other  parameter. 

As  an  example,  a  50;«  increase  in  rotor  inertia  makes  a  significant  change 
in  the  H-V  envelope  for  the  0H-6A  and  negligible  change  in  the  rating 
factors.  See  Table  A4  and  Figure  A-4(a).  On  the  other  hand  the  same  iner¬ 
tia  increase  for  the  AH-1G  yields  a  negligible  change  in  the  H-V  "nose" 
point  and  high  hover  he i yh t  out  a  significant  change  in  the  rating  fac¬ 
tors  occurs,  See  Table  4  and  Figure  A-4(b).  This  indicates  that  minimal 
benefits  can  be  attained  by  increasing  the  AH-1G  rotor  inertia  unless  a 
change  in  touchdown  body  attitude  is  permitted  so  as  to  utilize  the  bene¬ 
fits  of  the  inertia  change. 

One  might  argue  that  the  rating  factor  could  also  be  interpreted  as  a 
safety  margin.  This  safety  margin  could  serve  as  an  indicator  of  the 
margin  of  error  or  skill  level  for  the  pilot  in  conducting  the  auto- 
rotational  landing.  This  is  a  valid  argument,  and  the  significance 
of  increased  inertia  would  be  based  on  the  magnitude  of  the  rating  factor 

In  either  case,  the  combination  of  H-V  envelope  change  and  rating  factors 
provide  a  useful  guide  for  determining  the  effect  of  parametric  changes 
on  the  autorotative  capabilities  of  the  subject  helicopter. 

A  table  showing  the  effect  of  singular  parametric  changes  from  baseline 
configuration  of  the  0H-6A  and  AH-IG  is  provided  in  Table  A-4.  The  inter¬ 
pretation  of  this  table  would  be  as  follows. 
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Based  on  the  "optimum"  H-V  envelope  having  a  rating  factor  of  1.0,  the 
baseline  0H-6A  at  a  disc  loading  of  3.5  could  be  improved.  Review  of 
the  various  parameters  shows  that,  using  a  009.5  airfoil  section,  a 
5%  increase  in  the  low  hover  height  could  be  attained  with  slightly 
improved  error  margins.  This  airfoil  change  would  also  provide  a  14; 
reduction  in  the  "nose"  point  speed  and  the  high  hover  height  and  bring 
the  rating  closer  to  an  "optimum"  design  point.  Obviously,  other  design 
considerations  must  be  weighed  prior  to  making  such  a  change. 

Review  of  the  AH-1G  baseline  data  (Table  A-4)  shows  that  at  a  disc  loading 
of  5  no  one  parametric  change  could  significantly  improve  the  baseline 
H-V  envelope  but  rather  a  combination  of  parametric  changes  would  be  re- 
qui red. 

The  parametric  trends  shown  in  Table  A-4  are  only  valid  for  the  disc 
leadings  for  which  they  were  generated.  This  disc  loading  basically 
specifies  the  speed  range  for  rotorcraft  evaluation.  The  effect  of 
parametric  changes  in  other  speed  ranges  could  be  significantly  different. 
An  example  of  this  is  shown  in  Figure  A-5.  This  figure  shows  the  effect 
of  disc  loading  of  the  two  bas'eline  rotorcraft.  It  can  be  seen  from 
these  data  that  as  disc  loading  increases  the  H-V  envelope  boundaries 
increase.  The  sensitivity  of  airframe  drag  in  different  velocity  ranges 
could  be  significantly  different.  Other  major  parameters  are  similarily 
affected.  The  effect  of  increasing  altitude  has  a  similar  effect  on 
increasing  the  h-V  envelope  as  shown  in  Figure  A-G.  In  view  of  this,  any 
analysis  conducted  for  the  purpose  of  improving  or  designing  rotorcraft, 
for  minimizing  H-V  envelope  boundaries  must  be  conducted  for  the  intended 
operational  disc  loading  and  altitude  for  the  rotorcraft  in  question. 

A  parametric  trade-off  can  then  be  conducted,  for  the  design  point,  and 
a  listing  of  the  sensitivities  compiled  in  a  similar  manner  as  Table  A-4. 

The  parameter  sensitivity  list,  as  presented,  is  only  intended  for  use  as 
a  guide  of  design  considerations  for  improving  au uoroiati onal  character¬ 
istics  of  helicopters  and  is  only  applicable  for  the  disc  loadings  for 
which  it  was  developed. 
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HDP  CORRELATION  WITH  FLIGHT  TEST  DATA 


To  establish  the  confidence  level  of  the  Helicopter  Dynamic  Performance 
(HOP)  Program,  a  comparative  analysis  between  calculated  and  flight  test 
performance  was  conducted.  This  comparative  analysis  was  categorized  to 
cover  the  major  assumptions  used  in  the  analysis.  The  rotorcraft  used 
for  this  study  were  the  Hughes  0H-6A  and  the  Bell  AH-1G. 

Initially,  the  HDP  Program  was  used  to  evaluate  steady-state  hover  per¬ 
formance  for  the  subject  rotorcraft  at  various  wheel  heights.  The  results 
of  this  evaluation  are  shown  in  Figure  A-7.  The  fairings  used  for  the 
test  data  were  obtained  from  Reference  8  for  the  0H-6A  and  from  Reference 
9  for  the  AH-1G.  The  results  of  this  comparison  show  excellent  agreement 
between  theory  and  flight  test  at  all  wheel  heights  tested.  It  should  be 
noted  that  although  an  apparent  discrepancy  exists  for  the  0H-6A  at  high 
C  levels,  these  values  are  only  attainable  at  high  altitudes.  The  maxi¬ 
mum  power  levels  for  sea  level  standard  and  10,000  ft.  standard  are  shown 
on  this  curve  for  reference.  Based  on  the  power  limitations,  the  maximum 
deviation  between  calculated  and  flight  test  is  10  SHP.  This  is  regarded 
as  being  within  measuring  accuracy  and  test  scatter  in  the  basic  hover 
data . 

The  apparent  deviation  between  the  calculated  and  faired  curves  at  the  low 
Cw  values  for  the  AH-1G  can  be  accounted  for  by  referring  to  the  actual 
measured  test  data  used  in  establishing  the  comparison.  The  measured  oata 
for  the  low  Cw  range  is  included  on  Figure  A-7  for  reference.  Also  shown 
for  reference  ia  a  Cp  band  of  2  x  10~5  which  equates  to  55  SHP  at  sea 

level  standard  conditions. 

Program  verification  for  forward  speed  performance  was  conducted  with  tho 
HDP  program  set  in  a  steady-state  mode  for  the  speed  range  of  the  rotor¬ 
craft  in  question.  The  comparison  between  calculated  and  measured  test 
data  for  the  0H-6A  at  sea  level,  5000  and  10,000  ft.  density  altitudes  is 
shown  in  Figure  A-8.  The  test  data  used  for  this  comparison  were  obtained 
from  Figures  32,  33,  and  40  of  Reference  8. 

Comparisons  of  level  flight  performance  for  the  AH-1G  in  the  clean  and 
heavy  hog  configuration  are  shown  in  Figure  A-9.  The  test  data  for  these 
figures  were  obtained  from  Figures  9  and  4  (respectively)  of  Reference  10. 
The  calculated  level  flight  performance  shows  excellent  agreement  with  the 
measured  test  data,  thus  verifying  the  analytical  rotor  performance  model 

"ENGINEERING  FLIGHT  TEST  OF  THE  0H-6A  HELICOPTER  (CAYUSE)  PHASE  D,  Final 
Report,  USATEC0M  Project  No.  4-6-0250-C1,  USAASTA  Project  No.  65-37, 

April  1969. 

"ENGINEERING  FLIGHT  TEST  AH-1G  HELICOPTER  (HUEY  COBRA)  Phase  D,  Part  2 
Performance,  Final  Report,  USATEC0M  Project  No.  4-6-0500-01,  USAASTA 
Project  No.  66-06,  April  1970. 

’“ENGINEERING  FLIGHT  TEST  AH-lG  HELICOPTER  (HUEY  COBRA)  Phase  B,  Part  6 
Final  Report,  USATEC0M  Project  No.  4-6-0500-01,  USAASTA  Project  No.  66-06, 
November  1969. 
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in  level  flight.  It  is  interesting  to  note  that  a  7  sq.  ft.  aF  was  re¬ 
quired  above  the  AH-1G  clean  configuration  to  obtain  the  Heavy  Hog  con¬ 
figuration  performance.  This  is  in  agreement  with  the  flight  test  find¬ 
ings  as  reported  in  Reference  10. 

With  the  rotor  performance  model  verified  in  powered  flight,  a  further 
analysis  was  conducted  to  establish  its  validity  in  power-off  (autorcta- 
tive)  steady-state  flight.  The  results  of  these  calculations  are  shown 
in  Figure  A-10.  The  measured  test  data,  obtained  from  Figure  55, 

Reference  8,  for  the  0H-6A  and  Figure  15,  Reference  10,  for  the  AH-1G 
were  superimposed  on  the  calculated  curves. 

Review  of  the  0H-6A  comparison  indicates  a  high  degree  of  test  scatter 
in  the  data.  In  an  attempt  to  account  for  this  scatter,  several  computer 
runs  were  conducted  to  establish  the  sensitivity  of  the  descent  rate 
with  variations  in  trimmed  tip  path  plane  angle.  The  results  of  this 
study,  superimposed  on  this  figure,  show  a  100-FPM  variation  in  rate  of 
descent  per  1  degree  variation  in  tip  path  plane  angle  from  the  trimmed 
condition  (0H-6A).  In  view  of  this,  along  with  the  variation  in  the 
test  data  at  a  given  speed,  a  550-FPM  variation  at  35  knots  850-FPM 
variation  at  45  knots,  and  220-FPM  variation  at  55  knots,  it  was  judged 
that  the  program  output  was  well  within  the  measuring  accuracy  of  the 
test  data.  The  tests  conducted  with  the  AH-1G  demonstrated  a  higher 
degree  of  repeatability  and  the  calculated  results  are  within  ±1.5 
ft. /sec.  of  the  measured  data.  The  sensitivity  of  the  AII-1G  rate  of 
descent  with  tip  path  plane  angle  variations  from  trim  was  calculated  to 
be  100  FPM/0.75  deg.  A  tip  path  plane  angle  from  steady-state  trim. 

This  portion  of  the  comparative  analysis  covered  the  major  aspects  of 
the  rotor  performance  model.  The  agreement  shown  between  theory  and 
flight  test  barically  verified  the  assumptions  made  in  regard  to  ground 
effect,  number  of  blades,  blade  compressibility  and  stall,  and  the 
distribution  of  power  measurements  in  regard  to  induced,  profile,  parasite 
and  climb.  The  ability  to  represent  the  airfoil  data  as  a  cubic  equation 
in  a  mean  lift  coefficient/mean  drag  coefficient  relationship  was  also 
verified.  It  should  be  noted  that  the  OH-GA  has  a  four-bladed  rotor 
system  using  a  0015  airfoil  section  while  the  AH-1G  uses  a  two-bladed 
rotor  with  a  9.33  percent  symmetric  airfoil  section. 

To  verify  the  dynamic  portion  of  the  program,  three  autorotative  flight 
test  landings  conducted  with  the  AH-1G  were  simulated  with  the  HDP. 

These  landings  were  initiated  from  a  steady-state  flight  condition  with 
a  simulated  engine  malfunction  followed  by  a  1  to  2  sec.  collective  delay 
and  rul 1  autorotative  touchdown. 

The  first  condition  evaluated  was  a  power  cut  from  a  15-foot  hover, 
(representative  of  the  height-velocity  low  hover  height).  The  results 
of  the  calculated  flight  were  superimposed  on  the  test  data  traces  in 
Figure  A-ll(a).  The  calculated  results  show  excellent  agreement  with 


"HEIGHT  VELOCITY  TEST  AH- 1 G  HELICOPTER  AT  HEAVY  GROSS  WEIGHT  (LOW  ELEVA¬ 
TION),  final  Report,  USAASTA  Project  No.  74-19,  June  1974. 
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rotor  RPM  droop.  A  2-second  time  delay  was  used  for  this  flight.  The 
vertical  accelerations  and  displacements  along  with  the  collective  stick 
position  also  show  excellent  agreement  with  the  flight  test  test  time 
history.  The  flight  test  data  was  obtained  from  Figure  6  of  Reference  11. 
The  summary  sheet  of  the  completed  flight  path  is  shown  in  Tigure  A-ll(b). 

A  forward  flight  entry  condition,  72  knot  approach  at  195  feet,  (nominally 
representing  the  H-V  nose  point)  was  evaluated  and  the  results  are  shown  in 
Figure  A-12(a).  The  summary  sheet  for  the  calculated  flight  is  presented 
in  Figure  A-12A(b).  Evaluation  of  a  throttle  chop  flight  from  670  ft.  AGL 
(indicative  of  AH-1G  H-V  high  hover  height)  is  shown  in  Figure  A- 13 ( a )  with 
the  computed  flight  summary  sheet  shown  in  Figure  A- 13(b).  These  data  were 
obtained  from  Figures  3  and  1  of  Reference  10. 

Review  of  these  data  shows  excellent  agreement  between  theory  and  flight 
test  The  program  shows  some  variation,  however,  in  vertical  accelera¬ 
tion  and  airframe  pitch  attitude  during  the  mid-portion  of  the  maneuver 
i . e . ,  after  the  time  delay  and  prior  to  the  cyclic  flare.  The  program 
also  indicates  that  the  rotor  RPM  starts  to  drop  off  slightly  sooner  and 
at  a  slightly  higher  rate  during  the  final  collective  flare  than  the 
test  data  indicates.  The  reason  for  this  is  not  readily  apparent  as  the 
calculated  vertical  accelerations  and  airframe  pitch  attitudes  are  in 
nominal  agreement  with  measured  data. 

It  should  be  noted  that  the  true  airspeed  trace  shown  for  the  flight 
test  data  refers  to  flight  path  speed.  Although  some  discrepancy  exists 
between  calculated  and  measured  velocities,  excellent  agreement  between 
the  calculated  and  measured  horizontal  distances  was  obtained. 

Detailed  time  history  data  for  autorotative  landings  with  the  OH-6A 
rotorcraft  was  not  available  for  this  study.  In  view  of  this,  two  auto¬ 
rotative  landings  conducted  with  the  Sikorsky  CH-54B  were  used  for  this 
comparative  analysis.  While  these  flights  do  not  include  the  entry  into 
autorotation  they  do  pick  up  the  flight  during  the  steady-state  descent 
and  final  flare  maneuver.  Flights  at  a  nominal  sea  level  standard 
density  altitude  and  10,000  feet  density  altitude  were  used  for  this 
analysis.  These  data  are  shown  in  Figure  A-14. 

In  the  HDP  analysis,  C  „  refers  to  the  mean  blade  angle  required  over 
the  total  rotor  disc  for  generating  the  desired  thrust.  This  blade 
angle  is  referenced  to  the  tip  path  plane  axis  system.  1  he  relationship 
between  ‘yt,  and  the  collective  stick  position  requires  the  transformation 
from  the  tip  path  plane  axis  system  to  the  airframe  axis  system  along 
with  the  collective  rigging  of  the  rotorcraft  in  question.  The  HDP 
analysis  presently  does  not  incorporate  the  capability  to  conduct  this 
angle  transformation  or  provide  for  rotorcraft  rigging  inputs. 

However,  a  evaluation  of  the  transformation  equations  required  was  con¬ 
ducted  externally  from  the  program  for  an  axial  flight  condition  with 
the  AH-1G.  This  condition  was  selected  since  the  tip  path  plane  and 
body  attitude  were  constant  throughout  the  flight.  Based  on  the  rigging 
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assumption  that  full  collective  stick  travel  wi 1 1  produce  a  21-degree 
change  in  mean  blade  angle  for  a  neutral  CG,  th*i  transformation  of  075 
(average  blade  angle)  to  collective  stick  position  shows  excellent  agree¬ 
ment.  (Figure  A-ll).  In  conducting  this  evaluation,  the  HDP  program  was 
set  up  to  use  the  A075  input  option.  This  option  provides  the  capability 
of  specifying  the  mean  blade  angle  and  computing  .the  resulting  flight  path 
and  rotor  RPM  changes  with  time.  As  can  be  seen  in  Figure  A-ll,  the  time 
history  of  rotor  RPM  and  flight  path  shows  excellent  agreement  with  flight 
test,  thereby  verifying,  at  least  from  hover,  the  energy  concept  for  cal¬ 
culating  the  rotor  mean  blade  angle.  The  transformation  of  the  mean  blade 
angle  to  collective  stick  position  in  translational  flight  would  require 
the  details  of  the  control  rigging  relationships  for  the  rotorcraft  in 
question. 
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SPECIFICATION  COMPLIANCE 


The  0H-6A  and  AH-1G  autorotative  capabilities  with  regard  to  their  com¬ 
pliance  with  MIL  Spec.  MIL-H-8501A,  General  Requirements  for  Helicopter 
Flying  and  Ground  Handling  Qualities,  were  reviewed.  The  pertinent 
paragraphs  of  the  MIL  Spec,  for  this  evaluation  were  3. 5. 4. 4,  3.5.5. , 
3.5.5  1  and  3.5.7.  The  HOP  evaluation  for  this  specification  compliance 
check  was  established  using  HOP  time  history  output.  The  results  of  96 
autorotative  landings  are  summarized  in  Tables  A-l  and  A-2.  Presen¬ 
tation  of  the  detailed  time  histories  is  beyond  the  scope  of  this  effort. 

Paragraph  3.5. 4. 4  -  from  a  35-knot  autorotational  touchdown,  bring  the 
helicopter  to  a  stop  within  200  feet. 

HOP  calculations  indicate  that  this  criterion  cannot  be  met  with  a  skid 
type  gear  when  a  level  asphalt  surface  is  used. 

Paragraph  3.5.5  -  Enter  into  autorotation  at  all  speeds  from  hover  to 
maximum  forward  speed--Col 1 ecti ve  pitch  control  motion  has  been  delayed 
for  at  least  2  seconds  following  loss  of  power. 

HOP  analysis  indicates  that  a  2-second  time  delay  can  be  tolerated  for 
AH-1G  normal  operating  weights  up  to  75%  maximum  forward  speed.  This 
2-second  collective  delay  cannot  be  safely  tolerated  at  overload  oper¬ 
ating  weights  or  at  forward  speeds  in  excess  of  75%  maximum,  but  must  be 
reduced  to  1  second.  This  is  based  on  the  AH-1G  test  data  (Reference 
13)  and  the  HOP  analysis  for  these  data  shown  in  Figures  A-12  and  A-13. 
Again,  no  0H-6A  autorotation  time  history  data  was  available  at  the  time 
this  report  was  written. 

Paragraph  3. 5. 5. 1  -  Sudden  loss  of  power  with  collective  control  fixed, 
shall  not  produce  pitch,  roll,  or  yaw  attitude  changes  in  excess  of  10 
degrees  in  2  seconds,  except  that  at  speeds  below  that  for  best  climb,  a 
20-degree  yaw  in  2  seconds  will  be  accepted. 

HOP  analysis  indicates  that  both  aircraft  can  meet  the  intent  of  3. 5. 5.1 
for  the  pitch  and  roll  attitude  rates  within  the  limitations  of  permis¬ 
sible  collective  delay.  See  3.5.5  above.  The  HDP  program  does  not 
evaluate  yaw  attitudes. 

Paragraph  3.5.7  -  It  shall  be  possible--to  make  repeatedly  safe,  power-of 
autorotative  landings  (touchdown)  at  speeds  of  15  knots  or  less. 

HDP  analysis  indicates  that  neither  the  0H-6A  nor  the  AH-1G  can  comply 
with  this  specification  at  normal  operating  gross  weights. 


CONCLUSIONS 


This  study  was  conducted  to  establish  the  sensitivity  of  various  design 
parameters  on  the  autorotational  characteristics  of  helicopters.  The 
Helicopter  Dynamic  Performance  (HOP)  program  was  used  for  the  study. 
Based  on  the  results  of  this  effort,  it  is  concluded  that: 

The  sensitivity  cf  any  given  design  parameter  is  highly  dependent  on 
the  baseline  conditions  from  which  the  sensitivity  is  desired. 

The  sensitivity  of  a  given  parameter'  on  a  given  rotorcraft  cannot  be 
directly  used  as  typical  for  all  rotorcraft  even  when  the  same  baseline 
H-V  capability  exists. 

A  highly  generalized  effect  of  various  parameters  can  be  categorized,  as 
shown  in  Table  A-3.  In  using  this  table  the  designer  must  keep  in  mind 
that  these  parameters  are  highly  interdependent  and  the  changing  of  one 
might  not  be  effective  without  the  modification  of  others  which  would 
permit  attaining  the  desired  change. 

The  benefits  of  parametric  changes  on  autorotational  capabilities  must 
be  evaluated  on  an  aircraft  by  aircraft  basis. 

The  designer  must  also  consider  the  effects  these  parametric  changes 
might  have  on  the  overall  power-on  flight  envelope  for  the  intended 
design. 

The  sensitivity  of  parametric  changes  as  predicted  by  the  HDP  program  is 
based  on  control  input  requirements  which  maximize  the  aircraft  cap¬ 
abilities.  Pilot  capabilities  for  duplicating  this  maneuver  are  not 
considered.  The  resulting  calculated  H-V  envelope  must  then  be  con¬ 
sidered  as  the  minimum  possible. 

The  results  of  the  comparative  analysis  between  program  calculated  per¬ 
formance  and  flight  test  data  show  excellent  agreement  in  all  areas  in¬ 
vestigated.  In  view  of  this,  the  Helicopter  Dynamic  Performance  program 
can  be  used  with  a  high  degree  of  confidence  for  estimating  rotorcraft 
capabilties  in  power-on  and  emergency  flight  conditions. 
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RECOMMENDATIONS 


The  Helicopter  Dynamic  Performance  (HOP)  program  developed  under  this 
contract  provides  a  useful  tool  for  evaluating  rotorcraft  takeoff  and 
landing  procedures  and  capabilities.  Its  use  is  intended  to  be  a  guide 
for  evaluating  rotorcraft  designs  and/or  flight  procedures.  Its  use 
does  not,  however,  suggest  the  elimination  of  flight  tests  to  verify 
"real  life"  aircraft  and  pilot  workload  capabilities. 

It  is  recommended  that  this  program  be  used  prior  to  actual  flight 
testing  to  provide  some  insight  as  to  basic  aircraft  capabilities  and 
possibly  dangerous  conditions.  This  preliminary  effort  will  minimize 
the  time  and  cost  of  required  buildup  flight  test  as  well  as  minimize 
the  risk  in  conducting  them. 

Its  use  as  a  tool  for  evaluating  various  design  concepts  is  also  recom¬ 
mended.  Analysis  of  parametric  sensitivities  should  be  conducted  in  a 
manner  similar  to  that  used  herein.  The  change  in  H-V  envelope  from  the 
baseline,  variations  in  rating  factors ,  and  overall  effect  on  the  in¬ 
tended  flight  envelope  should  be  fully  evaluated  prior  to  committing  to 
final  design  or  configuration  changes. 

To  improve  the  versatility  of  this  program,  it  is  recommended  that 
additions  to  the  program  be  made  for  evaluation  of  wings  and/or  tail 
planes,  and  co-axial  rotors.  Further  development  for  coupling  the  HDP 
program  with  an  acoustic  analysis  to  provide  acoustic  foot  prints  or 
detectability  criteria  during  powered  takeoff,  landing  and  nap-of-the- 
eartn  maneuvers  should  be  considered. 

Further  analysis  should  be  conducted  to  better  establish  the  rela¬ 
tionship  of  the  mean  blade  angle  (  C^),  as  defined  for  the  energy  method, 
and  the  collective  stick  position  during  translational  flight. 

The  program  should  also  be  modified  to  permit  the  input  of  known  flight 
test  control  input  time  histories  for  correlation  with  flight  test  data. 

Additional  correlation  over  a  wider  range  of  helicopter  types  is  recom¬ 
mended.  This  correlation  could  lead  to  the  development  of  a  validated 
and  calibrated  autorotative  index  for  "rule-of-thumb"  evaluation  of 
helicopter  autorotative  capaDi 1 ities. 

The  HDP  program  can  also  be  used  for  effectively  optimizing  power-on 
takeoff  and  landing  procedures  for  small  field  operation  studies. 
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Figure  A-l.  Height  Velocity  Comparison  of  Baseline  Configuration 


(b)  EFFECT  OF  PARAMETRIC  VARIATIONS 
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Figure  A-8.  0H-6A  Level  Flight  Performance  (Clean  Configuration) 
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(b)  AH-  1G 


Figure  A-iO.  Amorotation  Performance 
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Figure  A- 12.  Time  History  of  Autorotational 
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(b)  Calculated  Data  (Forward  Speed  Approach) 
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figure  A-13.  Time  History  of  Autorotational  Landing. 


figure  A- 13.  Concluded. 


igure  A- 14. 


table  a-i.  parametric  variations  and  effects 


NOTE:  Run  No.  2  is  baseline  case. 


TABLE  A-2.  PARAMETRIC  VARIATIONS  AND  EFFECTS 


TABLE  A-3.  GENERALIZED  PARAMETER  SENSITIVITY 


Effect  on  H-V 

Parameter 

0H-6A 

AH-  1G 

Disc  Loading  (increased) 

degrades 

degrades 

Altitude  (increase) 

degrades 

degrades 

LG  (moving  aft) 

degrades  • 

degrades 

Rotor  Inertia  (increase) 

improves 

improves 

Allowable  Flare  Angle 
(decrease) 

degrades 

degrades 

Airfoil  (decrease  drag) 

i mproves 

improves 

Time  Delay  (increase) 

degrades 

degrades 

Initial  NR  (increase) 

improves 

improves 

Airframe  Drag  (increase) 

degrades 

degrades 

Airframe  Static  Stability 
(decrease) 

no  change  noted 

no  change  noted 

Vertical  Impact  Speed 
( i ncrease) 

improves 

improves 

Body  Altitude  at  Touchdown 
(decrease) 

degrades 

degrades 

Max.  Allowable  Cy/a 
(decrease) 

degrades 

degrades 

TABLE  A-4.  PARAMETRIC  VARIATIONS  FROM  BASELINE 
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10SE"  IN  SPEED 


LIST  OF  SYMBOLS 


a 

a' 

a" 


a 

TOT 

ANG 

b 

B 


C 


Rotor  lift  curve  slope,  acceleration 
Axial  interference  factor 

Radial  interference  factor 

Horizontal  Acceleration 

Rotor  disc  area 

Area  of  front  rotor  (tandems) 

Total  disc  area 

Total  projected  disc  area 

Tip  path  plane  angle 

Number  of  rotor  blades  on  rotor 

Rotor  efficiency 

Blade  chord 

Drag  coefficient 

Vertical  drag  coefficient 

Wetted  area  drag  coefficient 

Effective  blade  chord 

Centrifical  force 

Mean  lift  coefficient 

Available  power  coefficient 

Climb  power  coefficient 

Induced  power  coefficient 

Profile  power  coefficient 

Parasite  power  coefficient 

Stall  power  coefficient 

Torque  coefficient 


ft/sec 


ft/sec 

ft2 

ft2 


ftL 

ft2 

deg 


ft 


ft 

lb 
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Drag  torque  coefficient 
Thrust  coefficient 
Weight  coefficient 

Longitudinal  distance  between  tandem  rotors  ft 

Rotor  interference  factor 
Rotor  interference  factor  (hover) 

Rotor  interference  factor  (forward  rotor  on  rear  rotor) 
Rotor  interference  factor  (low  rotor  co-axial  system) 

Rotor  interference  factor  (level  flight) 

Rotor  interference  factor  (rear  rotor  on  front  rotor) 

Rotor  interference  factor  (upper  rotor  co-axial  system) 


Drag  lb 

Rotor  drag  component  in  climb  lb 

Profile  drag  lb 

Parasite  drag  lb 

Resultant  distance  main  rotor  centroid  to  airframe 
center  of  gravity  ft 

Distance  from  airframe  center  of  gravity  to  tail  rotor 
centroid  (Horizontal  for  pitch;  vertical  for  roll)  ft 

Vertical  drag  lb 

Blade  hinge  offset  ft 

Function  of 

2 

Equivalent  flat  plate  area  for  drag;  force  ft  ;  lb 

Unbalanced  force  in  longitudinal  direction  lb 

Unbalanced  force  in  lateral  direction  lb 

Unbalanced  force  in  vertical  direction  lb 

2  2 

Acceleration  of  gravity  -32.2  ft/sec  ft/sec 
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Gross  weigh;-. 


Initial  height 

H-V  .lose  point  height 

High  hover  height 

Lower  hover  height 

Horsepower  available 

Horsepower  to  climb 

Induced  horsepower 

Prof i le  horsepower 

Parasite  horsepower 

Stall  horsepower 

Wheel  height  at  arbitrary  point 
2  Mh 

Hub  moment  -  Cl  be  — 

Horsepower 

Horsepower  available  -  tip  rockets 
Horsepower  available 
Horsepower  to  accelerate  horizontally 
Horsepower  to  accelerate  vertically 
Horsepower  derived  from  rotor  RPM  change 
Horsepower  required 

Unsupported  power  loading 
Power  required  to  hover  OGE 

Airframe  inertia  about  longitudinal  axis 

Airframe  inertia  about  lateral  axis 

Airframe  inertia  about  vertical  axis 


ft-lb 


ft-lb-sec^ 


f t-1 b-sec 


ft-lb-sec 


Ut 


K 

KE 

M 

MB 


'l  ,90 


n 

N 

N 


P. 

l 


M 

QF 

R 


S 

S 

71 

SHP 

t 

T 


Rotor  inertia  about  rotation  axis 
Time  step  counter 
Kinetic  energy 
Mass 

First  flapping  moment  of  rotor  blade 
Advancing  tip  Mach  number 
Profile  drag  correction  factor 
Arbitrary  percentage  of  weight 
Aerodynamic  efficiency  in  climb 
Rotor  speed 
Power 

Airframe  pitch  acceleration 
Induced  power 
Airframe  pitching  mrment 
Torque  factor 

Rotor  blade  radius  (measured  from  center  of 
rotation) 

Airframe  rolling  moment 

Displacement 

Wetted  area 

Shaft  horsepower 

Time 

Thrust 

Thrust  front  rotor  (tandems) 

Thrust  main  rotor 

Thrust  rear  rotor  or  rocket  thrust 
Tali  rotor  thrust 
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ft-lb-sec^ 

ft-lb 

slugs 

slug-ft 


ft-lb/sec 

2 

deg/sec 

ft-lb 

ft 

ft-lb 

ft 


sec 

lb 

lb 

lb 

lb 

lb 


7  ~  T 


-  jtflillhtt*'..  .  <  ...  -  V.  ,  ,1,^*1 


L 


u 

Induced  velocity 

ft/sec 

u0 

Induced  velocity  in  hover 

ft/sec 

u 

Induced  velocity  ratio  u/uQ 

- 

V 

Veloci ty 

ft/sec 

VCR 

Critical  velocity 

ft/sec 

VH 

Horizontal  velocity 

ft/sec 

VR 

FI ight  path  velocity 

ft/sec 

vs 

Vertical  distance  between  tandem  rotors 

ft 

Vv 

Vertical  velocity 

ft/sec 

\ 

Vertical  velocity  design 

ft/sec 

VW 

Wind  velocity 

ft/sec 

vx 

Velocity  at  arbitrary  point 

ft/sec 

V 

Flight  path  velocity  ratio  V^/Uo 

- 

w 

Aircraft  weight 

lb 

Airframe  pitch  attitude 

Airframe  pitch  rate 

Airframe  pitch  acceleration 

Airframe  pitch  attitude  (previous  time  step) 

Rotor  overlap  angle  (tandems);  climb  angle 

Mean  profile  drag  coefficient 

Increment  of 

Wake  skew  angle  (tandems) 

Impressed  blade  pitch  angle 
Equivalent  linear  blade  twist 
Impressed  blade  pitch  at  0.75  R 
Rotor  inflow  ratio 
Ground  effect  parameter 
Rotor  tip  speed  ratio 
Horizontal  component  of  tip  speed  ratio 
Tip  speed  ratio  at  minimum  power  point 
Vertical  component  of  tip  speed  ratio 
Constant  3.14159.  .  . 

Ai r  density 

Air  density  at  sea  level  standard 
Rotor  disc  solidity 
Summation  of 

Rotor  induced  velocity  angle 
Angular  velocity  of  main  rotor 
Design  angular  velocity  of  main  rotor 


deg 

deg/sec 

2 

deg/sec 

deg 

deg;  deg 


deg 

rad 

deg 


1  b-sec'Vft4 
lb-sec2/ft4 


rad/sec 

rad/sec 
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LIST  OF  EQUATIONS 


This  section  lists  the  major  equations  used  in  the  Helicopter  Dynamic 
Performance  Program.  All  equations  are  numbered  as  referred  to  in  the 
text. 


hp  =  hp.  +,  hp  +  hp  +  hp  +  hp 


i  ro  rP 


2.  T  =  p*fr  (VR  +  u)  2u 


3. 


•  ],/2 


2pTrfr 


V  * 


-u  sin  (u  sin2aw-  u2  +  1/u2 ) (Reference  1) 


5.  u  =  [ - 1 - ]V2 

actual  rotor  2ttp(BR)2 


6. 


550  hp .  =  CT  3/2 


pirR2(nR)3  B 


1.  =  [A0/A]1/2 


8.  A/A  =  1  -  l^-s™Ycosy 

o  tt 


2  ,1/2 

-  i  LA-  1 

2  1 1  AT0T 


9-  pi  =2Pi  tr  ] 


10.  dfy  =  2  [ 


1 


1  _  Y  -  s i n yCQS y 

IT 


] 


1/2 


dfLF  -  dfFR  +  dfRF  =  2/0  +  PZcosZ0 


2 _ 2 


11. 


I'  .  .  =  ‘-5 T72cAu  [1  ♦  JL 

total  L 


‘(H) 


•2  CT/a  [r - 1 - , 

r  +  bm2  -  v 

3  r  97 


CP0  =  T~  0  +  nn2) 


hP„  =  °e6  ,„d2/,.Dx3 


o  (1  +  ng2)  ElB-ML 
8  ‘  550 


u--5±cJ  +  rV/2 


2pirRc 


u/u0  =  -  IT  i  0/4  (jj-)2  *  l]1 


*  «')B  ■  U(J  .  [ — 


2pnfr 


u(l  +  a')  =  [ - I _  ]l/2 


2dt t(BR  ) 


+  a')  =  550  hp.  /T 


24. 


[1.34Ct]  1/b 

F  =  1 - g — —  -  [0.143258e  +  0.035] 

«.  V/g  =  V/[  T_  7]V2 
0  2ttp(BR)^ 


26. 


ABu  =  xm(2/Ct)1/2  +  +  1]1/2  -  ♦  1] 

i  uT  tT 


1/2 


26  ( a  1 


aBu  =  0.0905u  [p-]1/2  +  +  1]1/2  -  [£^6974/  +  ] 

LT  ^LT  Ct 


1/2 


[1.34CT]1/b  .  1/2  2  , 

E  =  1 - g-1 -  +  0.0905m  [q-]  7  +  [fc~  +  1]  U 


27. 


-  [L|974m_  +  -j  -jl/2  _  £0  14325  0  +  0.035] 


28.  Mff  =  [VH2  +  Vv2]1/2/,iR 


29.  z1  0_(x-O.8) 


x3d  x  =  /  ^ 


8  X 


3. 

x  d  x 


30.  d  .  £|I 


31.  C 


F  3 


PP  2nR2 


32  hP0  ,  r  PnR2^R)3 
p  LPp  550 
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— ^ - .. TT,:.7T!I.\:.^a41.  M’C.'  x» 


/2Bwv  2  ^T__ 

33.  CT  cos  aTpp  cyj  +  C0v^U/,Uo  +  q^/2  2B2 

33(a).  <^-28*  (1  -^) 


34. 

hpc  s  (T 

cos  aTpp)Vv  (po/p)/550 

35. 

*CIIMB 

=  (CT  cos  aTpp)uv(p0/e) 

36. 

o 

ro 

II 

1  1  7 

CxJdx/  /  x  Jdr. 

0 

0 

37. 

_  b 

°e  rR 

/  1  Cx3dx/  f  1  x3dx  =  4b/*  /  1  C/R  x3dx 

0  0  0 

38. 

f(CQ/o. 

Ml,90*  6e)  *  [1750. (CD/o)2  -208.5  CD/a(3.6M2-3. 71 (M+l ) ) 

+(((CD/0  +  .01)2/. 000133)-. 752)(M-.7)+38.7M-21M2-9. 9] 
(1 .17+1 .21 75©e ) 


39.  m<-25 

f(CL/o-(CL/o)crit)=(l.-(.33+3.u)CL/o+.133u+.036)-(CL/o)crit 

40.  y> .25 

f(CL/a-(CL/o)crjt  )--CL/o-(Cyo)crit-(l  .08Cl/o-.071  ) 

41.  CL/o  -  (CL/o)LSL  <0.018 
f(CL/o-(CL/o)LSL)  =  3.0(CL/o-(CL/a)L5L) 
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42. 


CL/o-(CL/o)LSL>  0.018 


f(CL/o-(CL/a)L$L)  =  833.33(CL/o-(CL/o)LSL)' 

+  30.65(CL/o-(CL/o)LSL)  +  .336 

43. 

+  [f(CL/o-(CL/o)LSL)]a(CL/o)LSL), 

44.  hp  =  Cp  pitR2(qR)3/550 

KS 

45.  0tTPP ( K+l  )  ,  ANGq  +  (ANGmax-ANGq)  Sin  |>/2( 

V  -  V(K) 

a-j-pp(K+l)  =  ajpp  sin  [rr/2(-^-r- - )]  + 

Vz  x 

47.  SHP( K+l )  =  SHP(K)  +  HPMAK/( tf^yg.) 

48‘  Hpn  s  CV1100^3]2  '  r‘22)/at 


t(Kvl)-t 

49.  N  =  N  +  aNr  sin  0/2  (  — - £)] 

R0  R  tRATE 


50. 


a  - 


6  C 

(1+1 .5u2)— -  +  1 . 5 >.  (1-0. 5p2) 

Oa 

(1-p2+2.25u4) 


51.  a  =  -<J> 


u  v^CT 

52.  A  -  ~Uu s i na  -  -  —■ - 

1  * 


t(K+l)-t 

t - £)1 

RATE 

tan'^D/T) 


53.  <P  =  tan  *(y) 


54. 


r  _  0~u  +2 . 25m  )a-l .  5 A ( 1  *0 . 5ij^)  /oa\ 

t  (r> 


55. 


n2 *  I#i  -  —s—-^ 

M 


56.  F  =  TcosaTpp  sinaTpp  -D 
x  TPPL0N6  LAT 


57.  F  ■  Tsina^pp  C0SaTPP  -D 
y  p*L0NG  1KKifiT  y 


LAT 


58.  F  -  TcOSajpp  COSaypp  - ( GW+D. . ) 
2  PPL0NG  PPLAT  V 


59.  a  =  F/M 


6°.  VNEW  =  V0LD  +  adt 


61. 


AS  =  V0LDdt  +  1/2  adt2 


62  ■  ST0TAL  "£  o1  -s 

63'  W(X)  =  Vwcos(DIRx-HOG)  cos(DIRz) 

64-  W(Y)  =  Vwsin(DIRx-HOG)  cos(DIRz) 

65.  W(z)  =  Ywsin(0IRx) 

66  HPACC  =  AKE/550 
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I  ii  hi  1 1  mo  ■  mt 


1 


67.  h,  =  [ 


.  .11  i  n 

1.1  -V/V  ~]  "  J1 


68.  V  =  [ 


•  11V  (hx-  y,) 


(1-1V.-V) 


2  J  H 


69.  a^cos'1^, 


PM  s  ^TR^TR  ’  ^  ^R  TMR  +  s’n(-6^. 


71.  6  =  +  6t.x  -6q 


72.  I  =  I  +  w/gDc' 
WmR  ^CG  R 


73>  p*  =  Pm/1  =  B 
A  M  ^MR 


74-  Bt  =  6t-l  +  6tAt 


75.  6t  =  Btat  +  .5StAt  +  6t_1 


74.  (m)  6t  =  ( Bt_ +  0tAt)PDArtP 


76.  HPuns  =  HPjl(l-HPA„/HPjO 


AV 


77.  h 


W\ 


10  =  '  noo  hpuns  ^'2'24  v/ct/o) 


78.  C^/o  -  2/u|^j^  Cj/a 


2 

79.  „WN  -  0.95  [(CT/B)2M-  ]1/4 


80.  VCR  =  (2.809Vmin  +  5.613  CJj  -169.776)  ll- 


HP 


AV 


/  GW  VVn\ 


) 


81-  =  0.18Vcr2  +  199.0 


82. 


hhi 


•SyG- 

nR*" 


HP 


AV 


(' 


GW  VVq 
HP"  A+  '  550 


)  +  lll.v 


83.  hx  =hlQ  ♦  (0. 11/(1. 1-Vx/VCR)  -- H)(hCR  -h]o) 


84.  „x  -  0-(l-Vx/VCR)'/2)  (hM  -  hj 


hi  CR 
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